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We study the Coulomb blockade in a chaotic quantum dot connected to a lead by a single channel 
at nearly perfect transmission. We take into account quantum fluctuations of the dot charge and a 
finite level spacing for electron states within the dot. Mesoscopic fluctuations of thermodynamic and 
transport properties in the Coulomb blockade regime exist at any transmission coefficient. In con- 
trast to the previous theories, we show that by virtue of these mesoscopic fluctuations, the Coulomb 
blockade is not destroyed completely even at perfect transmission. The oscillatory dependence of 
all the observable characteristics on the gate voltage is preserved, its period is still defined by the 
charge of a single electron. However, phases of those oscillations are random; because of the ran- 
domness, the Coulomb blockade shows up not in the averages but in the correlation functions of 
the fluctuating observables (e.g., capacitance or tunneling conductance). This phenomenon may be 
called ^mesoscopic charge quantization" . 

PACS numbers: 73.23.Hk, 73.40.Gk 



I. INTRODUCTION 

The effect of Coulomb blockade§~i in chaotic quantum 
dots sets up an ideal stage for studying the interplay be- 
tween the quantum chaos and interaction phenomena in 
a many-electron system. By tuning the connection be- 
tween the leads and the quantum dot, one can study a 
rich variety of non-trivial effects. In the weak tunneling 
limit, discrete charging of the dot results in a sequence of 
sharp conductance peaks, which carry information about 
the chaotic motion of non-interacting electrons confined 
inside an almost closed doto. In the opposite limit of 
wide channels, charge quantization does not occur, and 
quantum chaos of free electrons in an open billiard may 
be studieoH. In a broad intermediate region, the charge 
quantization is gradually destroyed, and the chaotic elec- 
tron motion is affected by fluctuations of chacge of the 
cavity. The modern experimental techniquetrQ'Q allows 
one to continuously traverse between these regimes. 

The effect_of charging is conventionally described by 
HamiltonianH 

where C is the total capacitance of the dot, the dimen- 
sionless parameter Af is related to the gate voltage V g and 
gate capacitance C g by Af = V g /eC g , and Q is the dot 
charge. Usually, charging energy Ec is much larger than 
the one-electron mean level spacing of the dot, A. If the 
connection of the dot with the leads is weak and temper- 
ature T is small, T <C Ec, the charge is well quantized 
for almost all Af except narrow vicinities of the charge 
degeneracy points (half- integer Af) . The behavior of the 
differential capacitance of the cavity, d(Q)/dV g and of 
the conductance through the cavity is quite different for 
the system tuned to the immediate vicinity of charge de- 
generacy points (Coulomb blockade peaks), or away from 



those points (Coulomb blockade valleys). The statistics 
of the peaks can be related to tha,properties of a single 
electron energy and wave functioncl, so that the distribu- 
tion functions for these quantities can be extracted from 
the well known Random Matrix theory (RMT)Ba The 
transport in the valleys occurs by virtual transitions of 
an electron via excited states of the dotEil The statistics 
of the conductance in this case was recently obtained in 
Ref. |ll| and was confirmed experimentally in Ref. |l2|. 

All the aforementioned results were obtained neglect- 
ing quantum fluctuations of the charge of the cavity. 
These fluctuations grow with the increase of the coupling 
between the dot and leadtSl. Then, the difference between 
the peaks and valleys becomes less pronounced and even- 
tually instead of the peak structure, one observes only a 
weak periodic modulation^. Clearly, this modulation 
can be described neither by the properties of the single- 
electron wave function nor by the lowest order virtual 
transitions via the excited states. 

The case of almost perfect transmission of a one- 
channel point contact connecting the quantum dot 
with the leadj-see Fig. |l|, was analyzed by Matveevtil 
and Flensbergjij in the framework of an effective one- 
dimensional Hamiltonian. Employing the bosonization 
technique, they showed that the Coulomb blockade dis- 
appears completely if the transmission coefficient of the 
point contact is exactly unity r — and A = 0. The ex- 
plicit dependence of the differential capacitance of the 
system on Af, and on the small reflection coefficient 
< |r| 2 <C 1 was obtained by MatveevtJ. It is important 
to emphasize that the Coulomb blockade in this situation 
is non-perturbative in charging energy effect and it can 
not be revealed in the standard Hartree-Fock or Random 
Phase approximations. 

The properties of a quantum dot connected to reser- 
voir by a channel was analyzed|-4p a series of papers 
of Biittiker and collaboratorai^O. They have used 
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the Random Phase approximation to calculate the fre- 
quency dependence of the linear response of the current 
I through the channel to bias Ve luJt applied to the reser- 
voir: / = G{lo)V . In this approach, the admittance 
C = — Im(dG/du>)\ u —>Q, coincides by construction with 
the thermodynamic capacitance of the non-interacting 
electrons. The quantum corrections and mesoscopic fluc- 
tuations of this quantity can be then analyzed by using 
the distribution of_the Wigner delay times of the non- 
interacting systemO. This approach is perfectly valid 
for a large number of channels, but for a single channel, 
there is no any parameter justifying it. We will see below 
that the results obtained by well controlled procedure are 
significantly different, see Sees. [v|, and [v| 




FIG. 1. Schematic view of a quantum dot connected to a 
lead. Periodic orbit "A" encounters the entrance to the dot 
once, riA = 1, and periodic orbit "B" does not encounter the 
entrance, tib = 0. 

In this paper, we account for both the strong quan- 
tum charge fluctuations, and the chaotic electron mo- 
tion within the dot. We will show that backscattering 
of electrons from the walls of the dot into the channel 
connecting it to the lead results in residual Coulomb 
blockade oscillations of observables with the gate volt- 
age. In the limit of perfect channel transmission, the 
relative magnitude of the differential capacitance oscilla- 
tions is x/A/Ec and (A/E c ) In 2 (E c /A) for the spinless 
and for spin 1/2 cases respectively. If the second lead is 
attached to the cavity by a weak tunnel junction with 
conductance Go <C e /h, the two-terminal conductance 
G can be measured. The average value of the conduc- 
tance is suppressed by the Coulomb blockade, but it re- 
mains finite even at zero temperature, (G) ~ GoA/Ec- 
Fluctuations of the conductance are of the order of its 
average. This resembles the result for the elastic cotun- 
neling in weak coupling regimePH 11 !. However, the depen- 
dence of the conductance of the gate voltage is no longer 
a sequence of deep valleys and sharp peaks, but rather a 
weakly oscillating function. 



For a finite reflection coefficient in the channel (|r| 2 ^ 
0) , we found a new contribution, in addition to the aver- 
aged differential capacitance calculated in Ref. [l4]. This 
contribution is fluctuating and provides, in particular, 
the dependence of the differential capacitance on the 
magnetic field. 

The paper is organized as follows. In Sec. [l]|, we qual- 
itatively discuss the mesoscopic fluctuations of the dif- 
ferential capacitance for spinless electrons. Section III 
is devoted to the formulation of the model and deriva- 
tion of the effective action representation. We will also 
discuss the conditions of applicability of the model. Sec- 
tion IV describes the bosonization procedure. Calcula- 
tion of the ground state energy and differential capaci- 
tance is performed in Sec. [v]. The tunneling conductance 
in a strongly asymmetric setup (one channel is reflection- 
less, and the other junction is of conductance Go <C e 2 /h) 
is studied in Sec. [V^ . Our findings are summarized in 
Conclusion. 



II. QUALITATIVE DISCUSSION 

Let us consider first a completely opened channel (r = 
0). In the limit A — > 0, the electron charge of the dot 
varies with the gate voltage as (Q) t^-A^V, to assure the 
mini mum of the electrostatic energytj. The interaction 
fll.lD depends only on the number of electrons crossing 
the dot-channel boundary. Therefore, the properties of 
the ground state can be characterized by the asymptotic 
behavior of the wave- functions far from the entrance to 
the dot. This behavior is described by the scattering 
phase, and at low energies can be understood from the 
following qualitative argument. 

Entrance of an additional electron with energy e (all 
the energies will be measured from the Fermi level) into 
the dot requires energy Ec- Therefore, the electron may 
spend in the dot time of the order of Ti/ Ec, and then the 
extra charge of the dot has to relax. There are two pro- 
cesses, that lead to the relaxation of the charge: i) elastic 
process where the same electron leaves the dot; and ii) in- 
elastic process where some other electron is emitted from 
the dot. At low energies the probability of the inelastic 
process is small as (e/Ec) 2 , by virtue of the smallness 
of the phase volume. (The last statement assumes the 
Fermi liquid behavior at low energies and, as we will see 
later, is valid only for spinless one-channel case.) There- 
fore, we may consider only elastic process. The same 
consideration is applicable also to an electron leaving the 
dot. Thus we conclude, that the low energy properties 
of the system can be mapped onto the dot effectively de- 
coupled from the channel, and the phase of the scattering 
amplitude from the entrance of the dot is given by the 
Friedel sum rule 



S = 7r(Q)/e = ttN. 



(2.1) 



Equation (2.1) can be applied to electrons incident from 
inside the dot, as well as to electrons incident from 
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the channel. The outlined description resembles closely 
the Nozieres description-of the unitary limit in the one- 
channel Kondo problemliS. 

The outlined above qualitative picture based on the 
introduction of scattering phase S is somewhat intuitive; 
it will be verified by a calculation in Sec. Here in- 
stead of rigorous proof, we demonstrate that this scheme 
reproduces the result 



E g (Af) ~ \r\E c cos 2irAf 



(2.2) 



obtained by MatveevQ for the ground state energy 
Eg (AT) of spinless electrons in the limit of zero level spac- 
ing in the dot. Then, we apply the scheme to find the 
corrections to the ground state energy arising from a fi- 
nite A. Those corrections will result in the mesoscopic 
fluctuations of the ground state energy. 

We start with considering the limit A = 0. First, we 
put also |r| = and calculate the density of electrons in 
the channel p(x). Then we take into account the scatter- 
ing potential V{x) that generates r ^ 0, in the first order 
of perturbation theory, 



E g {N) = J dx P {x)V{x). 



(2.3) 



As we discussed, the Coulomb interaction leads to the 
perfect reflection of electron at low energies; wave- 
functions have the form ipk(x) = cos (k\x\ — 5), with 
phase shift 8 given by Eq. (2.f). It leads to the 
Friedel oscillation of the electron density p(x) = 
J2 VF \k-k F ]<E c I^M^)! 2 , where v F and k F are the Fermi 
velocity and Fermi wavevector respectively. We obtain 



p{x) 



E c 

— cos(2fc F H - 25), \x\ < vp/Ec; 
v F , 



sin(2fc F |x| - 2(5) 



|x| > v F /E c 



(2.4) 



Here we omitted the irrelevant constant part of the elec- 
tron density. Substituting Eq. (2.4) into Eq. (2.3), as- 
suming the magnitude of the potential around x = 
smaller than v F /Ec, and using the standard expression 
|r| = \V(2k F )\/v F , we obtain formula (|J). Here V(k) 
is the Fourier transform of the potential V(x). 

Having verified the suggested scheme for the caseA = 
0, we proceed with evaluation of the ground state en- 
ergy of a finite dot connected to a reservoir by a per- 
fect channel. According to the above discussion, the 
channel is effectively decoupled from the dot due to 
the charging effect even though r = 0. Therefore, we 
have to relate the ground state en ergy of a closed dot 
to the scattering phase 5 of Eq. (2.1). For a chaotic 
dot, this problem is equivalent to finding a variation of 
the eigenenergies by introduction of impurity potential 
V(r) = (\/ttv) S(r) tan 6, where v is the averaged density 
of states per unit area. The relevant contribution to the 
ground state energy is given by 



E„ 



E 

-E o <5«<0 



(2.5) 



where £j are the eigenenergies measured from the Fermi 
level p. As soon as the scattering phase changes by w, 
one more level enters under the Fermi level. Evolution 
of the energy levels with changing 5 is shown schemat- 
ically in Fig. ||. The position of the level £i(<5) satisfies 
the gluing condition 



&(<5 + 7r)=& +1 (<5). 



(2.6) 



From Eqs. (2.5) and Eq. ( |2.6| ) we see that the ground 
state energy depends almost periodically on 5: 



Eg(6)=Eg(6 + ir)+0(A). 



(2.7) 



As we will see below, the amplitude of the oscillation of 
the ground state energy with 6 is much larger than the 
mea n lev el spacing and we will neglect the last term in 
Eq. A 



In order to estimate the amplitude of the oscillations, 
we recall that the correlation function of the level veloc- 
ities is given byc3 



(ds <H 



A 

n 



(dseidse 3 



(2.8) 



where /3 = 1,2 for the orthogonal and unitary ensem- 
bles respectively, a nd ( . . .) stands for the ensemble av- 
eraging. Formula ( |2.8| ) can be easily understood from 
the first order perturbation theory. At <5 <C 1, we have 
Ci(S) « ej(0) + (5/^)1^(0)]. For a chaotic system the 
exact wave functions ipi can be presented in the form 



itMo)r 



A 



where the area of the dot A appears due to the normal- 
ization condition, and bi characterizes the fluctuations 
of the chaotic wavefunctions. In accordance with the 
Porter-Thomas distributionO, 



(bi 
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FIG. 2. Evolution of the energy levels of the quantum dot 
with the scattering phase 5. 

Estimating the mesoscopic fluctuation s of the ground 
state energy ( |2.5| ) with the help of Eq. (2.8), we obtain 
for 5 < 1 



{[E g (§)-E g (0)Y 



(2.9) 



AM- 

7T 



E 



-Eo<ti,t}<0 



AE C S 2 . 



As we have already explained, energy (2.5) is a periodic 
func tion of S with period n. On the other hand, for 8 < 1, 
Eq. (|1^) is valid. Therefore the characteristic amplitude 
of the oscillations is of the order of ^/EqA, and it is plau- 
sible to assume that the correlation function of energies 
at two different parameters Afi , M 2 takes the form 

(E g (M) E g (A/2)) « A£ c cos 2vr (A/i - A/2) , (2.10) 



where we use Eq. (2.1). It is important to notice that 
the variation of the energy of the ground state is much 
larger than the mean level spacing A. This observation 
enabled us not to consider in particular the variation of 
the chemical potential with changing 6, because it would 
generate a correction of the order of the level spacing 
only. 

In order to expla in the functional form of the correla- 
tion function ( 2.10 ) and make our argumentation more 
precise, it is instructive to evaluate the shift of the_energy 
levels starting from the Gutzwiller trace formulaEj. The 
energy of the ground state is given by 



E n 



deN(e)K(e/E c ), 



(2.11) 



where e is measured from the Fermi level, and N(e) is the 
integrated density of states, N(e) — J2i — e 0- Here 



K(x) is some function which decays rapidly at x > 1, so 
only the states which can be described within the Fermi 
liquid theory contribute into the 5-dependent part of the 
ground state energy [cf. Eq. ([[1])]. 

The integrated density of states caxL.be expressed as a 
sum over the classical periodic orbitsa: 



N(e,5) = Re^i?j(e)exp 



Me) 



2irij5 



(2.12) 



Here Rj is the weight associated with j-th orbit, Sj is 
the reduced action for this orbit. The last term in the 



exponent in Eq. ( 2.12 ) characterizes the reflection from 
the entrance of the cavity, and rij is the number of such 
reflections for j-th orbit, see Fig. "J. We have omitted the 
mean val ue of iV(e) which is independent of the phase 



shift (|2.lD. 

Integrated density of states ( 2.12| ) is a strongly oscil- 
lating function of energy which vanishes upon ensemble 
(or energy) averaging. However, it contributes to the 
fluctuations of the density of states: 



{N(e 1 ,5 1 )N(e 2 ,S 2 ))=2ReJ2\Ro 



(2.13) 



exp 



-(5 J (e 1 )-^( e2 )) + 2m i (5 1 -J 2 ) 



and we neglected the energy dependence of the pre- 
exponential factors Rj because it occurs on the energy 
scale of the order of the Fermi energy. In the double sum 
over the periodic orbits, arising in (2.13), one can retain 
only diagonal termsEil because different orbits have dif- 
ferent actions; the non-diagonal terms oscillate str ongly 
and vanish upon the averaging. Factor of 2 in Eq. fl2.13| ) 
originates from the fact that, the electron trajectory, j, 
and the trajectory time reversed to j have the same ac- 
tion. 

In order to calculate correlation function of meso- 



scopic fl uctua tions of ground state energies ( 2.11 ), we 
use Eq. (2.13) and expand action as Sj(e±) — Sjfa) = 
(ei — £2)tj, with tj being the period of j-th orbit. After 
integration over energies ei ; 2, we find 

(EgiSjEgfa)) = (2.14) 

2El £ |^ | 2 A (^p-) cos [2 nj (6 1 - 5 2 )} 



V ft 



is a function decaying 



where K (y) = f_ dxK(x)e lxv 
at y > 1. 

Coefficients \Rj\ 2 have a very simple physical mean- 
ing, and arc related to the classical probability P(t)dt 
to find a periodic orbit with a period within the interval 
[t;t + dt]: 



(2.15) 
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In the same fashion, we obtain from Eq. ( [2.14 ) 



(2.16) 



(¥) ^(t)cos[2n(«W 2 )] . 



The classical probability P n (t) differs from P(t) by sat- 
isfying an additional constraint: the corresponding peri- 
odic orbits are reflected from the dot entrance exactly n 
times; £ n P n (t) =P(t) . 

In general, P n (t) depends on a particular shape of the 
dot. However, if the motion is chaotic, P n (t) becomes 
universal, 



Pn(t) 



(tA/h) n tA/h 

1 ' 



(2.17) 



for the periods t much larger than the ergodic time 
Ti/Et- Energy Et associated with the time scale at 
which the classical dynamics becomes ergodic is the coun- 
terpart of the Thoulcss energy for the diffusive system. 
Typically Et ^ Ec, therefore we adopt approximation 
E T ^> E'c > A. According to Eq. (|1J), the character- 
istic period of the semiclassical trajectory is of the order 
of %/ Ec, which is much smaller than the Heisenberg time 
h/A. Therefore, when calculating the oscillatory part of 
the correlation function of the ground state energies, it 
suffices to keep in Eq. (2.16) only the contribution of 



the trajectories reaching the entrance only once. As the 
r esult , we obtain the expression similar to our estimate 

( B ), 



(E g (Ai) Eg (A/2)) = aAE c cos 2tt (A/i - JV 2 ) , (2.18) 



where we have used Eq. (2.1). The numerical coefficient 
a = 4j f °° dxK(x) depends on the particular form of 
function K, and cannot be found within the simple con- 
sideration. 

Let us now discuss the correlation of the ground state 
energies as a function of magnetic field. The magnetic 
flux threading a periodic orbit adds a phase <j> — j ; = 
AjH/^Q to th e action in each term j in the Gutzwiller 
formula ( |2.12 ), where Aj is the directed area under the 
trajectory, H is the applied magnetic field, and $0 is the 
flux quantum. Correspondingly, formula (2.14) is modi- 
fied to 



(Eg(S 1: H 1 )E g (S 2 ,H 2 ))=2E c J2\Rj\ 2 K 



E c tj 



E x A 
$0 



' \ <, )S f E^k. ) ( . as [2„ ( (A, - A,)] . (2.1!)) 



In analogy with Eq. (2.16), we transform Eq. (2.19) to 
(E g (6x,Ex)E g (5 2 , E 2 )) = -£V / dA -P n (t; A) x 

K (^) cos (M) C os (^) co S [2n(5x-S 2 )} , (2.20) 



$0 



where P n (t; A)dA differs from probability P n (t) by one 
more constraint: the directed area swept by a trajectory 
lies within the interval [A; A + dA] . 

In a chaotic system, the probability P n (t; A) factorizes: 

P n (t;A)=P n (t)p(t;A). (2.21) 

Here P n (t) is defined by Eq. (|2.17|) , and the distribution 
function of the areas is Gaussian: 

1 



p(t;A) 



E T t 



exp 



A 2 



A(A 2 {t)) 



(A 2 (t)) = -j-A 2 . 



(2.22) 



The formula for (A 2 (t)) shows that the typical area under 
the electron trajectory differs from the area of the dot A, 
and grows as \ft. This law is applicable at the time scale 
exceeding the ergodic time %/Et, and reflects the time 
dependence of r.m.s. of the random winding number for 
the trajectory of an electron bouncing off the walls of the 
dot. 

As we already discussed, the characteristic time an 
electron spends in the dot is fr/Ec- The characteristic 
area accumulated during this time is y/ErjEc- A mag- 
netic field produces an appreciable effect if a flux pene- 
trating through this area is of the order of 4> . Thus, the 
correlation magnetic field is controlled by the charging 
energy: 



E c 



$0 



E c 



(2.23) 



A V 2ttE t 

Using Eqs. (|Jo|)-(||J|) , we find the correlation function: 



E 2 



{E g (Afx,Hx)E g (M 2 ,H 2 )) = E C A £ A E [^j cos2™, 



7=+.- 



(2.24) 



where we introduced the short hand notation E± = 
Ex ±E 2 and n = Ai — Af 2 - Calculation of the exact form 
of dimensionless function Ae{x) = L dy e~ xy K(y), and 
of the numerical coefficient in Eq. ( |2.23| ) requires more 
involved treatment which is a subject of the following 
sections. 

Equations ( [2. IS ) and (2.24) constitute the main quali- 
tative result of this section. We were able to demonstrate 
the oscillations of the ground state energy with the ap- 
plied gate voltage. The phase of those oscillations is ran- 
dom, so that the oscillations can be revealed only in the 
correlation functions. Unfortunately, these simple quali- 
tative arguments arc not sufficient for finding the precise 
form of the correlation functions. Moreover, the assump- 
tion of the Fermi liquid behavior is valid only for the spin- 
less electrons. It is known that the loWifinergy behavior 
of the s = 1/2 electrons is equivalents to that of the 
two-channel Kondo problem in its strong coupling fixed 
point displaying a non-Fermi liquid behavior. Quantita- 
tive study of the system in this case will be presented 
later, see Sec: VB. 
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III. THE MODEL 

The main difficulty of the problem is in the non- 
perturbative nature of the Coulomb blockade effect. 
Derivation of an effective one-dimensional model is our 
first ste p in overcoming this difficulty. The interaction 
energy (1.1) depends only on the total number of elec- 
trons in the dot. The change of this number is asso- 
ciated with electrons propagating through the channel. 
Because the dynamics of the channel is one-dimensional, 
the charging effects of the system can be considered on 
the basis of a one-dimensional HamiltoniariiirEJ. How- 
ever, original problem was at least two-dimensional, so 
backscattering of the electrons by the walls of the dot 
can not be accounted for by the one-dimensional Hamil- 
tonian. Instead, of an effective Hamiltonian, we were 
able to find an effective action which depends only on 
the electron variables of the one-dimensional channel. If 
there were no interaction, such an approach would have 
no advantage; however, in the presence of interaction it 
becomes very powerful. The interaction will be ex actly 
accounted for by means of bosonization, see Section [TV]. 

Electrons are backscattered into the channel by the 
walls of the dot at rand om ti mes, therefore the action 
we derive in Subsection III A has a non-local 
random term 



m time, 
can be treated per- 



This term, however 
turbatively by virtue of the small parameter, A/E c <C 1. 
With the help of the action, calculation of the correlation 
functions of energies and differential capacitances can be 
performed by the standard diagrammatic methodsLl At 
energies less than Ep the averages become universal. In 
this regime, it is also possible to formulate thejnodel 
starting fr om th e Random Matrix HamiltonianE3, see 



Subsection IIIB 



T he a pplicability of description of the interaction by 
Eq. (1.1), i.e., of the const ant-interaction model, is dis- 
cussed in Subsection III C . We will show that the cor- 



rections to this description are of the order of l/g, where 
g ^> 1 is the dimensionless conductance of the dot. 



A. "Conventional" formulation. 

We start with the Hamiltonian of the system, 

H = H F + H C , (3.1) 

where Hp is the Hamiltonian of non-interacting elec- 
trons, 



H F = dr 



— W f Vi/> + (-fi + U(r)) ipty 
2m 



(3.2) 



The potential U (r) describes the confinement of electrons 
to the dot and channel. 



The interaction Hamiltonian He is given by Eq. (1.1) 
and the charge of the dot is 



9. 

e 



dot 



drip'ip, 



(3.3) 



where the integration is performed within the dot. Of 
course, the boundary separating the dot from the lead is 
not defined. However, this ambiguity can be absorbed 
into the definition of dimensionless gate voltage Af. 

For the purpose of the evaluation it is more convenient, 
however, to change the definition of the charge. Notic- 
ing that the total number of particles in the system is 
an integer number which can be added to the parameter 
AT without affecting any periodic in Af observables, we 
write 



channel 



drip* if) 



(3.4) 



To calculate the ground state energy, we start with the 
thermodynamic potential, 



In (ire"^) , 



(3.5) 



where temperature T = 1//3. 

We evaluate the trace in two steps, Tr . . . = TriTr 2 . . ., 
where 1 and 2 indicate the fermionic operators belonging 
to the channel and dot respectively. Because a ll th e in- 
teraction is attributed to the channel, see Eq. (3.4), the 
charging energy operator is not affected by the summa- 
tion in the dot, and can be omitted in the intermediate 



formulas. The non- interacting Hamiltonian (3.2) can be 
presented as Hp = Hi + H% + H12, where H\ and H 2 
are the noninteracting Hamiltonians of the channel and 
of the dot respectively, and Hyi connects the dot with 
the channel. Thus, we write 



Tr 2e -PH F = T -i3(h 1+ h 2 +h 12 ) = 



(3.6) 



g-jSQag-ZSHiy e 5 Jgdn J^<1t 2 {H 12 (t 1 )H 12 (t 2 )) 2 

Here H 12 (t) = e T ^ 1+fl2 ) H 12 er T ^ 1+fl2 ) is the interac- 
tion representation of the Hamiltonian connecting the dot 
and lead, T T stands for the chronological ordering, and 
f2 2 = — TlnT^e - ^ 2 is the thermodynamic potential of 
non-interacting electrons in the dot. Averaging (. . 
over the electronic variables of the dot is defined by the 

relation: (. . .) 2 = e /3r22 Tr 2 ^e _/3J ^ 2 . . . Thermodynamic 

potential Q 2 does not depend on Af, and it will be omit- 
ted. 



The operator Eq. (3J3) depends only on the electron 
variables o f th e channel. The evaluation of the last fac- 
tor in Eq. (3.6) is performed in Appendix |a[ This yields 



-(H 12 (t 1 )H 12 (t 2 )) 2 = -^(n;0)V(r 2 ;0) x 



dydy'<j)(y)(j>{y')dl xl Q(Ti - r 2 ;r,r' 



(3.7) 



Am 2 
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where ip(r;x) = e THl ip(x)e~ THl are the one dimensional 
fermionic operators of the channel in the interaction rep- 
resentation, i>(r) = t), and Q is the exact Matsub- 
ara Green function of the closed dot subjected to the zero 
boundary condition. The wave function 4>{y) describes 
the transverse motion in the single-mode channel, and 
the coordinates x, x' in the derivative of the Green func- 
tion Q are set to +0. 

The detailed behavior of the Green function Q de- 
pends on the particular shape of the dot. It is conve- 
nient to separate the fluctuating part of the Green func- 
tion, Q = Q + Q, and to combine the proportional to Q 
sample- independent par t of {Hvi{ti)Hii{t2))2 with the 
Hamiltonian H\ in Eq. (3.6). 

We rewrite Eq. (3.7) in the form 



{H 12 {t 1 )H 12 {t 2 )) 2 = (H 12 {t 1 )H 12 (t 2 )) 2 - 

2^(t i; 0)V>(t 2 ;0)L(ti-T2), (3. 

where kernel L is given by 

1 



L(r) 



Am 2 



dydy'mWKx'Q^T,*')- (3-9) 



One can check by a direct calculation that 

oc Ti^e-^ 10 , 



e -PH 1T e \ f Q <lTl J^dT 2 {Hl2(Tl)Hl2(T 2 ))2 



(3.10) 

where Tr> stands for the trace of the one-dimensional 
fermionic operators on the positive half-axis, and 



dx 



— W> f VV> - fiib^ib 
2m 



(3.11) 



is the one dimensional Hamiltonian defined on the whole 
real axis. The proportionality coefficient in Eq. ( 3.10| ) 
does not contribute to any observable quantity and we 
omit it. We substitute Eq. (3.8) into Eq. ( |3.6|) , use 
Eq. ( |3.10| ), restore the charging energy, see Eq. (|3.4D , and 
obtain 



Tre"^ oc Tr e^ ff °T T e 



-s 



(3.12) 



One-dimensional effective Hamiltonian is given by 



Hn — H 



ID 



E c Af 



tp^tp : dx 



(3.13) 



where : . . . : stan ds fo r the normal ordering. The effective 
action S in Eq. ( |3.12| ), 



S = 



d Tl dr 2 L (n - r 2 ) ^(n; 0)^(r 2 ; 0), (3.14) 



has kernel L defined by Eq. (3.9). If the electrons with 
spin are considered, the summation over spin indices is 
implied in the above formulas. 



As we are interested in the dynamics of the system at 
energies much smaller than the Fermi energy, we can lin- 
earize the spectrum of one dimensional fermions. Writing 



ip(x) = 



ipL(x) + e 1 i>n(x), where i/'l and ipn are 



the left- and ri ght-m oving ferm ions respectively, we ob- 
tain from Eqs. ( |3.1l|) and (|3.13|) : 



Hq = ivF / dx 



{^Ml - ^ R d x ^ R ) (3.15) 



wher e vp is the Fermi velocity in the channel. In 
Eq. ( |3. 15 ) the fermionic fields are assumed to be smooth 
on the scale of the Fermi wavelength. The action has the 
following form in terms of the left- and right-movers: 



S = 



dridr 2 L {t\ — t 2 ) x 



(3.16) 



[V'i(n) + <Mn)] [iPl{t 2 ) + tPr(t 2 )] . 

Finite reflection in the channel can be taken into account 



by adding to the Hamiltonian (3.15) one more term: 

cH bs = \r\v F (4(0)^,(0) + ^(0)Vl(0)) , (3.17) 

where \r\ 2 <C 1 is the reflection coefficient. The thermo- 
dynamic potential Q (J\f) can be found from 



Q(Af) = -TlnTr (VK^+^TVe^ 



(3.18) 



The differential capacitance Cdiff (A/ - ) of the system is 
then given by 



c m (AO = c l 



l d 2 n 

E^W 2 



(3.19) 



Equations ( |3.15| ) and ( |3.17| ) were first suggested in 
Refs. 

So far we succeeded in reducing the original problem 
to the effective one-dimensional problem, where all the 
features of the chaotic motion of electrons in the dot are 
incorporated into the non-local in time action. The ac- 
tion fluctuates strongly from sample to sample, and we 
should study the statistics of these fluctuations. 



1. Statistics of L(t). 

It is convenient to use the Lchmann representation for 
the function L(r) of Eq. (Oh: 



sinh [irT(t + it)] ' 

(3.20) 
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where t he r etarded and advanced kernels L R,A are given 
by Eq. (3.S), with Q replaced by the exact advanced and 
retarded Green functions Q R ' A respectively. It is well 
known that the averaged products of the type Q R Q R and 

Q A Q A vanish, and the products of the retarded and ad- 
vanced Green functions can be expressed in terms of the 
classical propagators - diffusons V D and Cooperons V , 



g R i (t 1 ;T+,r+)g A 2 (t 2 - 1 r-,r-) = 27nsS(t 1 +t 2 ) x (3.21a) 

T{v+, r-)f(r+ r-)V D (h; R X ,R 2 ) ; 

0* (ti; r+ r+)G A 2 (h;r^,r 2 ) = 2nu5(h + h) x (3.21b) 

T{v+^)T{v+,t 2 )V c (ti;Ri, R 2 ) , 

J 7 (ri,r 2 )=Im— = J (k F \r 1 - r 2 |)-J \ k F \ri - 7tr 2 |J , 

where Jq(x) is the zeroth Bessel function, = Ri±rj/2, 
and k F is the Fermi wavevector. Point IZr is the coor- 
dinate of the image charge crea ted b y the charge in the 
point r, so that the propagators ( 3.21 ) satisfy proper zero 
boundary conditions. Green functions here are ta ken a t 
different values of magnetic fields H\,H 2 . In Eqs. (3.21), 
v = m/2ir is the density of states per unit area. Equa- 
tions (3.21) are valid if the arguments of the Green func- 
tions are close to each other pairwise: |ri )2 | must be much 
smaller than the elastic mean free path for a diffusive 
dot, and much smaller than the dot size for a ballistic 
dot. Boundary is assumed to be smooth on the scale of 
the Fermi wavelength. 

Retarded classical propagators in the diffusive dot sat- 
isfy diffusion-like equationscJ 




V + — 

c 



Ai - A 2 
Ai + A 2 




= (5(Ri-R 2 ), 
(3.22) 



where D is the diffusion coefficient, and the vector po- 
tentials Ai >2 are defi ned so that V x A li2 = ifi )2 . For 
a ballistic dot, Eqs. ( 3.22j ) should be substituted by the 
corresponding Liouville (or, t o be more precise, Perron- 
FrobeniuSj) equation, and Eq. (3.21) should be somewhat 
changedEJ. However, in the universal limit considered in 
this paper, there is no difference between the ballistic and 
diffusive dots. 

The universal limit corresponds to a large time scale 
at which the semiclassical electron orbit covers all the 
available phase space. At such a time scale, the classi- 
cal probabilities no longer depend on the coordinate and 
acquire the form: 



V 



D.C 



(3.23) 



where 6(x) is the step function, A is the area of the dot, 
and the decay times associated with the magnetic field 
are given by 



D.C 



= E 7 



' H 



$1 T $2 



(3.24) 



Here $o = e/ch is the flux quantum, $i )2 = AH\ 2 are 
the fluxes through the dot corresponding to the fields H\ 
and H 2l and the Thouless energy is of the order of 
fiD/A for a diffusive dot, and of the order of hv F /\fA 
for a ballistic dot. 

The correlation functions of the retarded and advanced 
parts of the kernel L can be expressed, with the help 
of Eqs. ( |3.9| ) and (3.21), in terms of the diffuson and 
Cooperon. The kernel L, see Eq. (3.9), depends on the 
Green function at coinciding spatial arguments. There- 
fore, both pairings leading to the diffuson and Cooperon 
upon averaging, should be take n into account. In the 
universal regime, see E q. (f3~23 ), integrals in the trans- 
verse direction in Eq. ( |3.9[) can be calculated using the 
normalization condition J dy(f) 2 (y) = 1. As the result, we 
find 



• i 2 )0(tl 



(3.25) 



where A = l/(uA) is the mean level spacing of the dot. 
Averages of the type L R L R vanish. Because we are in- 
terested in times smaller than the Heisenberg time h/A, 
the higher moments can be decoupled by using the Wick 
theo rem a nd the pair correlation functions are defined by 
Eq. (HI). 



B. Random matrix formulation. 



We start from dividing the entire system into two 
parts, the leads and the dot. In general, the Hamilto- 
nian H of the system can be represented as 



H = Hr 



Hr 



H, 



LD- 



The Hamiltonian of the leads is of the form 



3' 



(3.26) 



(3.27) 



where we linearized the electron spectrum in the leads, 
and measured all the energies from the Fermi level. In- 
dex k is the longitudinal momentum in a mode propa- 
gating along the channel connected to the dot, and index 
j = 1, . . . , N labels these modes (summation over spin 
indices is implied whenever it is necessary). For the sake 
of simplicity, we assume the same Fermi velocity in all the 
modes, the general case can be reduced to Eq. (3.27) by 
the corresponding rescaling. Hamiltonian Ho = H n +Hc 
of the dot consists of the non-interacting part 
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a,/3 



(3.28) 



and the inter actio n term He, which is described by the 
Hamiltonian (1.1) with the charge Q given by 



Q 



(3.29) 



(In this subsection, we reserve Greek and Latin letters 
for labeling the fermionic states in the dot and in the 
leads respectively.) For definiteness, we restrict the dis- 
cussion to the case of orthogonal ensemble, generaliza- 
tion to ot her cases is straightforward. Elements Ti. a ,/3 in 
Eq. ( 3.28| ) form a real random Hermitian matrix TL of size 
M x M , (M — ► oo), belonging to the Gaussian ensemble 



P(TL) oc exp 



4A 2 M 



TrH' 



(3.30) 



where A is the mean level spacing near t he ce nter of the 
band. Finally, Hamiltonian Hld in Eq. ( 3.26 ) describes 
the coupling of the dot to the leads, and has the form 



Hld = J2 ( w aj*l>y>k<i + h-c.) , 



(3.31) 



where the coupling constants W a ,j is a real MxN matrix 
W. 

Let us list the needed averaged q uanti ties correspond- 
ing to the Gaussian ensemble Eq. ( 3.30 ). The averaged 
density of states is given by the semicircle law 



p(e) = Tr<S (e-H) = Re 



We will need only properties of the system at ener- 
gies e much smaller than the width of the band AN/tt, 
and we will neglect the energy dependence of the aver- 
aged density of states. Average of the Green functions 
g R ' A (e) = (e-H± M)^ 1 have the form 



n R,A . r 



MA' 



(3.33) 



Random matrix counterparts of the diffuson and 
Cooperon propagators (3.21) can be written as 



2tt 



+ V c {uo)8 ai 5 f 3s) , 
1 1 



MA 



M-iu + 0' 



(3.34) 



where Q = Q — Q. Averages of the type Q R Q R and 



Q A Q A vanish. Formulas (3.32) - (3.34) have the accuracy 
~ 1/M and they neglect the oscillatory w-dependences 
on the scale of the order of A. This accuracy, however, 



is sufficient for us because, as we already discussed, the 
relevant results are contributed by the energy strip of the 
width E G > A. 

Our purpose now is to deriv e the effective action theory 
similar to Eqs. (3.15) - (3. IS) starting from the random 
matrix model. Before doi ng so , let us review some useful 
properties of the system (3.26) in the absence of the in- 
teraction, Ec = 0. In this case, electron transmission at 
energy e is completely characterized by N x N scattering 
matrix S'(e): 

5(e) = 1 - 2mvWi [e-H + invWW^} ~* W, (3.35) 

where v = 1/(2itvf) is the one-dimensional density of 
states in the leads. Coupling matrix W can be repre- 
sented in the form 



W = 



AM 



-UOW, 



(3.36) 



where U is an orthogonal M x M matrix, W is a real 
N x N matrix and O is M x N matrix, O a j — Saj^ 1 < 
a < N. Because the distribution function ( |3.26| ) is 
invar iant under rotations TL — > UTLW , ma trix U in 
Eq. (13.36 ) can be omitted. Substituting Eq. (3.36) into 
Eq. ( [3.3E ) a nd pe rforming ensemble averaging with the 
help of Eq. (3.33), we obtain for the average scattering 
matrix: 



S 



1 - W^W 

i + w^w' 



(3.37) 



At W^ W = 1 the average scattering matrix vanishes, 



see Eq. (3.37). It indicates that matrix S belongs to the 
circular ensemble (corresponding to the regime of "ideal 
contacts"). Deviation of W^W from a unit matrix can 
be attributed to the scattering on the contacts between 
the leads and the dot. This scattering is described by 
unitary symmetric 2N x 2N matrix 



Sr.= 



with N x N matrices r, r 
1 - W^W 



ft 
r' 

t being defined by 



(3.38) 



1 + W^W' 

t = W 



r' = -W- 



1 + 



1 + w^w 



(3.39) 



The explicit relations between the coupling matrices W 
and the scattering matrix in the contacts (3.38) were first 
obtained by BrouweJ23. 

Now we turn to the derivation of the effective action. 



For technical reasons, [see di scussi on above Eq. (3.4)], we 
replace the charge operator ( 3.29 ) with 



(3.40) 
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After this replacement, the Hamiltonian of the system 
becomes quadratic in fermionic operators of the dot, so 
that this part of the system can be integrated out: 



Tine 



-PH 



(3.41) 



e -P(HL+Ho) e -p£l n j, e 3 J dT 1 dTv(H LD {T 1 )& LD (T2))r> 

Here Hld{ t ) is the interaction representation of the cou- 
pling operator Hld and averaging over Hamiltonian of 



the dot is defined as 



D 



Tr D (e-^ 



Thermodynamic potential flo = TlnTre~' 3 - ff ™ of the 
closed dot is independent of the gate voltage AT and it 
will be omitted. 

The average in the last facto r in Eq. (|3.4l|) is calcu- 
lated with the help of Eq. (3.31) and of the definition of 
the Matsubara Green function for the closed dot: 



Q a p(r) = -<T t ^ o (t)^(0)> d = £\ 



1 



iu> n — TL 



Ct0 



(3.42) 



where i/> q (t) 



ip a e' 



-H„t 



, ip( T ) = V^(~ r )j and 



u> n = ttT(2ti + 1) is the fermionic Matsubara frequency. 
The result is 



(Hld(ti)Hld{t 2 )) 2 



(3.43) 

^ i>kuh(n) [W^g(n - r 2 )w] jij2 V^ 2j2 (r 2 ). 

We separate the aver aged part of th e Green function 
Q = Q + Q, use Eqs. ( 3.33 ) and ( 3.36 ) , and obtain the 
Fourier transform of the kernel in Eq. ( ft.44| ), 



wg(iu n )w^ = 



; sgnw„ - 



XM -w^g{mo n )w. 



(3.44) 

The first term in the RHS of Eq. ( |3.44 ) does not contain 
any information about the dot, and its frequency depen- 
dence is the same as of the Green function of free chiral 
fermions. It is, therefore, possible (and very convenient) 
to transform this part of the action to the Hamiltonian 
form by introducin g fict itious fermionic fields bk.j, j = 
1, ...TV. Then, Eq. ( 3.41 ) acquires the form 



Trne-^ oc Tr b (e"^'"^^ 



(3.45) 



the omitted Af- independent proportionality coefficient is 
irrelevant. The effective Hamiltonian -ff e // is given by 



3,k 



— E K,ii M i'i»^A + 



h.c. 



+ 



ki,k 2 ;ji,j 2 



E c 



(3.46) 



where Wj 1 j 2 are the elements of the Hermitian matrix w 
defined as 



w 



1/2 



(3.47) 



Action S has the form 

r 

ki,k 2 ; 



5 = 4 V / dTidT2lpk 1 ,h(n)L jd2 (Ti-T 2 )'tpk 2 ,j 2 (T2), 
fc, t„. ,-. Jo 



(3.48) 



where the kernel L is a N x N matrix given by 
AM - 



L(r) 



(3.49) 



Equation ( 3.45 ) can be easily checked by tracing out 
fermions b and using the relation 

E ( T -r b kuji{ T )h 2 ,h( Q )) = i^S Jlj2 J2e ZbJnT sgnuj n . 

k\,k 2 



Hamiltonian H e f t can be rewritten in a more familiar 
form. Introducing the Fourier transform of the fermionic 
fields ipj(x) = y]i. e ~ ikx j >i.k and bj(x) = J2k e ~ ikx bj,k, 
we obtain from Eq. (3.46|) 



H, 



eff = >"F 



E 



dx 



{fid^i + b]i ),b,) 



+ 



h.c. 



+ 



(3.50) 



We are interested in the case of almost open ideal con- 
tacts, || w — 1 || <C 1. In this case, it is natural to change 
the variables and reveal the small parameter of the per- 
turbation theory. Introducing left- and right-moving 
fermions, 

i>j{z) = i J L,j(x)6(-x) + ip BJ (-x)0(x), 
bj(x) = i [ip R j(-x)9(-x) - tp Ltj (x)6(x)] 



(the ambiguity of this definition at the origin should be 
resol ved a s "0(0) = [VK+0) + ip{— 0)] /2), we obtain from 

Eq. (pbl) 
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H eff =H Q + H bs ; (3.51a) 

/oo 
dx H)\ ,()., i/i.j - tp^ jd x ij) R A + 

/(• f £ J dx: tPI^lj + ^r^r.j ■ A'j ; (3.51b) 



(3.51c) 



S = J2 d Tl dT 2 L l3 (n - r 2 ) x (3.51d) 

[i>L,i{n] 0) + $fl,i(Ti;0)] [Vfcjfa; 0) + Vfl,j(r 2 ;0)] . 

Here we neglected the terms related to the discontinu- 
ities of the fermionic field at the origin which induces 
higher order terms in r, and approximated reflection ma- 
trix r i=a 1 — w, as follows from Eqs. (3.47) and (3.39) for 
|| r ||<C 1. Within the same approximation we can put 
W = 1 in Eq. (|3.49| ). Formulas (3.51) are analogous to 
Eqs. ( 3.15 ) - (|3.18D derived in a previous subsection. 



Fin ally, we have to study the statistics of the kernel 
( |3~49| ) which can be expressed in terms of its adv ance d 
and retarded counterparts by Lehmann f ormu la (3.20). 



Performing time Fourier transform of Eq. (3.34), and us- 
ing v — 1/(2ttvf), we obtain 



(LM 1 )L^ s (t 2 )) = ^-S(t 1 +t 2 )e(t 1 )x 



2tt 
{SisS 



(3.52) 



Si r 5j s } . 



For the one-mode lead this result agrees with the zero 
magnetic field version of Eq. ( 3.25| ). 

The statistics of the kernel at different magnetic fields 
can be obtained by adding purely imagi nary Hetmitian 
matrix to the original matrix Ti. in Eq. ( 3.28| )BEll This 
would lead to the result analogous to the exponential de- 
cay in Eq. (3.25). We will not describe details of such 
calculation here, and refer the rea der , t o the extensive 
literature on parametric correlationsESca 



C. Applicability of the model. 

So far, we w ere using a very simple model of the in- 
teraction (1.1) which ascribes all the interaction effects 
to the variation of the number of particles in the dot. 
However, a natural question arises: what is the accuracy 
of this approximation? One may even think, that the ef- 
fects considered in this paper are completely washed out 
by remaining interaction terms which we neglected. 

The pur pose of this subsection is to show that the sim- 
ple model (1.1) of interactions in the dot leaves out only 
small, oc 1/g, effects, where g = Et/A is the dimen- 
sionless conductance of the cavity (Et is the Thouless 
energy). For a diffusive dot in the metallic regime, and 
for a ballistic non-integrable dot the conductance is large, 



g ^s> 1 . Mesoscopic charge quantizati on i s adequately de- 
scribed by the model of interaction (LI), as long as the 
j> 1, and the number of modes propagating to the lead 
is much smaller than g. In other words, ergodic time of 
the dot Ti/Et should be much larger than the escape time 
of electron from the dot. 

The electrons in the dot are described by the Hamilto- 
nian H = H n + Hi„ t , where the non -inter acting part of 
the Hamiltonian, H n , is given by Eq. (3.25). The validity 
of the random matrix theory for H n for the energy scale 
smaller than Thouless energy was proven in Ref. [2l] for 
chaotic systems and in Ref. |3(] for diffusive systems. 

The general form of the interaction Hamiltonian is 



Hi, 



2 ^ 



(3.53) 



In this subsection, we will write explicitly the spin indices 
a for the fermionic operators. The interaction Hamilto- 
nian (El) corresponds to the approximation of the ma- 
trix V by 



V « EcSatSp. 



7" 



(3.54) 



Our goal now is to show that all the other matrix ele- 
ments as well as the mesoscopic fluctuations of matrix 
elements (3.54) are small. Some-aL-these calculations 
already appeared in the literatureErE3, however, we will 
present brief derivation to make this paper self-contained. 

The easiest way to study the statistics of the one elec- 
tron wave- function </> Q (r) is to relate them to the Green 
function and then use Eq. (3.21). By definition of the 
retarded and advanced Green functions we have 

Gf (n , r 2 ) - Gf (n , r 2 ) = 2m ^ <\> a (n )<j> a (r 2 )<5(e - e a ), 

a 

(3.55) 

where we assumed no magnetic field for simplicity. At 
given e nergy e only one function contributes into the sum 
in Eq. ( 3.55| ), so that the statistics of the Green functions 
is related to that of the wave- functions. Furthermore, it 
is known that there is no correlation between level statis- 
tics and wave function in the lowest order in 1/g, see e.g. 
Ref. 



30 



so we can neg lect t he level correlations and aver- 
age (^-function in Eq. ( 3.55 ) independently. As the result, 
we can estimate 

Q (ri)0 Q (r 2 ) » A [G^(n,r 2 ) - G« (r ls r 2 )] . (3.56) 
2m 



Now we can use Eq. ( 3.55| ) and ( 3.2l| ) to study the aver- 
age of different momenta of the matrix elements 

Vap^s = J dr 1 dr 2 V(r 1 - v 2 )(j) a {r 1 )(j)p{T 1 )4> r y{T2)(j)s (r 2 ). 

Aver aging this matrix element itself with the help of 
Eq. (3.56), we obtain 
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Vc 



v. 



(0) 



+ v r 



(3.57) 



First term in Eq. (3.57) comes from the product of the 
averaged Green functions 



ImG A (ri,r 2 ) 



ik(ri-r 2 ) 



FS, 



and it is given by 

,(0) 



(3.58) 



V a p~/8 = Ec8 a s8^ + FA (SaySps + SaP^s) , 

Ec = -jp J dY 1 dr 2 V(r 1 - r 2 ), 
F = v(V(k)) FS . 



Here A is the area of the dot, (. . .) fs stands for the av- 
eraging over directions of the wave vector on the Fermi 
surface, v is the averaged density of states per unit area 
and per one spin in the dot, and V(k) is the Fourier 
transform of the two-particle interaction V(r). 

Charging energy E c in Eq. ( 3.58| ) is related to the 
zero mode of the interaction potential. This mode can 
not redistribute the electrons within the dot and that 
is why it is not screened (the redistribution of the elec- 
trons between the dot and the leads is taken into ac- 
count by the model). As the result, Eq is much larger 
than the mean level spacing. On the other hand, coeffi- 
cient F includes only non-zero modes which are perfectly 
screened, V{k) = V Q (k)/[l + 2uV (k)], where V (k) is the 
bare potential. (The use of the static screening here is 
possible because the screening is established during the 
characteristic time of the plasmon propagation through 
the dot, which is much smaller than Tl/Et-) Therefore, 
we estimate F < 1/2, so that the last two terms in- 
troduce a correction only of the order of level spacing, 
and may be neglected. (We will not consider here the 
case of th e attractive interaction when the third term in 
Eq. ( |3.58 ) renormalizes to infinity due to the interaction 
in the Cooper channel.) 

The second term in Eq. ( |3.57 ) originates from the prod- 
uct of the retarded and advanced Green functions, see 



Eqs. ( 3.21 ). In the absence of magnetic field, diffuson 
and Cooperon propagators coincide, and their spectral 
expansion for a diffusive system is 



1 



£ f;iri)Mr2 \ (3.59) 



7k #0 



where 7^ and / M (r) are the corresponding eigenvalues 
and eigenfunctions. For a diffusive system, 7^ = DQ^, 
where D is the diffusion constant, and wavevectors 
depend on the shape of the system. For a rectangu- 
lar diffusive dot of the size L x x L y , one finds Q 2 = 
ir 2 (n 2 /L 2 + n 2 /L 2 ) with n x , n y > being integer num- 
bers. For chaotic systems, 7^ are the eigenval ues o f the 
Perron-Frobenius operator. Zero mode in Eq. (3.59) cor- 
responds to the conservation of the number particles, and 
all the other modes describe the relaxation of any ini- 
tial inhomogeneous distribution function by the virtue of 



classical chaotic dynamics, Rej^ > 0. If the system is 
integrable, or there are some additional symmetries of 
the system, other zero modes appear, however, we dis- 
regard such a possibility and consider only diffusive or 
classically chaotic systems. 



Using Eqs. (3.21), (3.5£) and (3.59), and taking into ac- 
count that all the energies are smaller than the Thouless 
energy (or, in other words, the lowest non-zero eigenvalue 
of diffusion or Perron-Frobenius operator), \e a \ -C 71, we 
find 



V 



(1/9) 



[2Fi5 a s5pj + (F 2 + Fx) (SaySps + S a p5js)] , 

(3.60) 



where the dimensionless conductance of the system is de- 
fined as 



.9 = Re^, 



(3.61) 



and is assumed to be muc h lar ger than unity. Dimen- 
sionless coefficients in Eq. ( 3.6C ) are given by 



Ft 



5521 H SL, ,= 1,2; 



7T ' — ' 7,, 



(3.62) 



F? = u{V(k)) FS 

F 2=»f rfridr 2 V (n - r 2 ) /;(n)/ M (r 2 ). 

For the screened interaction potential, coefficients F^ 2 = 
1/2 and, therefore, iq. 2 are of the order of unity for 
chaotic systems and of the order of (l/Air 2 ) hx{L/l) for 
the diffusive dot. Here L is the size of the dot and I 
is the transport elastic mean free path. Thus, we have 
shown that the corrections to the average matrix element 
Q3.54 ) are parametrically small for the metallic regime. 

Now, we wish to show that the fluctuations of the 
matr ix el ements are small. Indeed, with the help of 
Eq. ( 3.21 ), we obtain for a generic (i.e., with no pair- 



wise equal indices) matrix element: 



A 



(YaPiS) — C ( — 



(3.63) 



The numerical coefficient c for the diffusive system is 
given by 



1 



(£) E{2(^) 2 + (^) 2 } Re 



7 M #0 



7,- 



[where coefficients F^ 2 are defined in Eq. ( 3.6 2| ) ] and it 
is of the order of unity, so that the matrix elements are 
small at g 3> 1. For chaotic systems, the expression for c 
is more cumbersome, but still have a similar stru cture. In 
the case of "diagonal" [in the sense of Eq. ( 3. 60) )] matrix 
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elements, the average in Eq. ( 3.63| ) should be replaced by 
the corresponding variance. 

The ma in conc lusio n of this section is contained in 
Eqs. ( 3.6C| ) and (3.63). These equations cle arly show 



that the Coulomb blockade type interaction (|l.l|) is a 
parametrically justified description for the dynamics of 
the system at energies smaller than the Thouless energy. 

Closing this subsection, let us m entio n numerical works 
that have been performed recentlyOa. These papers ad- 
dressed either dirty diffusive systems with small number 
of electrons!^, g ~ 1, or classically localized statesEJ and 
are not relevant for the metallic regime j > 1 we are 
dealing with. We believe that the large quantum dots 
studied in Refs. belong to the metallic regime. 



IV. BOSONIZATION PROCEDURE 



Equations (3.15) -(3.18) reduce the initial system con- 
sisting of a dot and a single-mode channel to the effective 
one-dimensional model. T o trea t interaction in the model 
(the second term in Eq. ( 3.15 )), we follow Ref. U, 

and 

use the bosonization tec hniqu e. In the bosonic variables, 
the entire Hamiltonian ( 3.1 5| ) becomes quadratic. The 
price for this convenience, is a strongly non-linear form 
that the backscattering terms acquire (in the language of 
left- and right movers, those are the terms oc ip^ipR in 



Eqs. (3.16) and ( 3.17] )). Fortunately, the typical value of 
the kernel (3.9) is small (~ A/Ec), and this enables us 
to use the perturbation theory, which will be presented in 
Sees. |v| and VI. In this section we present bosonization 
procedure in a form most suitable for our purposes. 



2vr 



(4.3) 



With the help of Eq. (gj), Hamiltonian ( j3~l~5| ) can be 
bosonized as 



Ec 
8^ 



dx 



( d<fii 



V dx 



dip 



[^l(0) + ^?(0) + 27tAA] 



dx 

2 



(4.4) 



The relations (|4.2a|) and (|4.2b|) ensure the fcrmionic 
commutation relations wit hin the species tpL and ipR. 
The commutation relation ( [4.2c|) guarantees the anticom- 
mutation relation of ipL with ip R and the commutation 
relations 



^U/jR + W L i>L<lx; tp RL (y) 



A L (yn-y)- 



Backscattering Hamiltonian (3.17) takes the form 
\r\v F 



H bs = 



ttA 



•cos [0l(O) + $n(0)] 



(4.5) 



Bosonized version of the effective action Eq. (3.16) is 



i r 

S=7TT d Tl dT 2 L ( n - t 2 ) j?(Ti)f}(r 2 ) x (4.6) 



2ttA 



A. Spinless electrons 



One-dimensional fermionic fields can be presented in 
the formEj 



^r{x) 



V 



n 



itp L {x) 



(4.1) 



where bosonic operators are taken at the origin x = 0. 
Majorana fermion rj does not enter the effective Hamil- 
tonian, and therefore it is not a dynamical field. Its role 
in the effective action is to take care of the difference in 
the definition of the operation of chronological ordering 
for the fermionic and bosonic operator. Equality 



(T t j7(ti)77(t 2 )) = sgn(n - r 2 ), 



(4.7) 



where A is the high energy cut-off of the order of the 
Fermi wavelength, and i) = if , t) 2 = 1 is the Majorana 
fermion, its significance will be discussed later. One- 
dimensional bosonic fields (pl,r(x) satisfy the following 
commutation relations: 



[<fiL{x),ifiL{y)] = -in sgn(x- y); 
[<Pr(x), ipniy)} = in sgn(x - y)\ 
[ip R (x),(f L (y)] = -in. 



(4.2a) 
(4.2b) 
(4.2c) 



It is e asy to check, using Eqs. (4.2), that the fermionic 
fields (4.1) obey the standard commutation relations. 
The expressions for the densities of left and right movers 
are 



and Wick's theorem preserves the defini tion of chrono- 
logical ordering for fermions in Eq. (3. IS). 

It is convenient to separate the part of the bosonic 
sector non affected by the Coulomb interaction and in- 
troduce new field tp+, tf—, $ with the commutation re- 
lations 



[<p+(x),<p+(y)] = -in sgn(x - y); 
YP-ix),^-^)] = -in sgn(x - y); 



[ip + (x),tp_(y)] = 0; 



(4.8a) 
(4.8b) 

(4.8c) 
(4.8d) 



We express operators (Ej) in terms of new fields (fhq) as 
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(pL{x) = t= 7rA/; 



<Pr(x) = 



V2 
V2 



(4.9) 



where c -number nAf is incorporated into the definition of 
t he f ield. It is easy to see that the commutation relations 
(4.2) are preserved. In the new variables, the Hamilto- 
nian (4.4) is independent on the gate voltage (i.e., on 
AO: 



( d<f A 



\ dx 



( d<P- 



\ dx 



(4.10) 



All the A/"-dependence is transferred now to the backscat- 
tering term in the Hamiltonian: 



£t rVF 
H bs = — — cos 

7TA 



V2(p+(0) - 2ttN 



(4.11) 



and to the action S which is contributed by the return 
trajectories of electrons after multiple scattering within 
the dot: 



1 f p 

S = — r / d,T 1 d,T 2 L (n - T 2 ) v{n)fi(T2) X 
ZwA J 



(4.12) 



exp 



. y-(ri) - (p-(T 2 ) 
V2 



exp 



.$(n) - $(r 2 ) 

1 7= 

V2 



7= \- -, — 7rA/ 

\/2 4 



^+(72) 
V2 



-ttJV" 



If one neglects such trajectories altogetherEl, then the 
Coulomb blockade oscillations apparently vanish in the 
limit r = (no backscattering in the channel). However, 
if return trajectories (i.e., a finite value of A/Ec) are 
taken into account, the Coulomb blockade oscillations ex- 
ist, even if r = 0. The A /~-dependence of the action ( 4.12| ) 
and Hamiltonian ( 4.11 ) clearly shows that the period of 
the oscillations does not depend on the details of the sys- 
tem. This periodic dependence is a direct consequence of 
the discreteness of the electron charge. 

We are going to develop perturbation theory in S and 
Hbs ■ Every order of the perturbation theory is expressed 
in terms of the correlators of t he b osonic field governed 
by the quadratic Hamiltonian (4.10). The necessary cor- 
relation functions are 



23_ (t)=<T t 0_(t)0_(O)); 
23 + (r) = (T t + (t)0 + (O)); 

23 $ (r) = (T T *(r)*(0)); 

23 $+ (t) = (tMt)0+(O)), 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



where averages are calculated with respect to Hamilto- 
nian Hq and all the bosonic fields are taken at x = 0. 

Standard calculation presented in Appendix |^ yields 
for T <C Ec and r > X/v F ~ l/e F : 



23$ (r) - 23$ (0) = - 



A 



dxe x In 



nT 



sin ( + t ) ttT 



T ] 7rT 



sinh' 



/ 2tt 2 Tx \ 



vf I sin7rTr| 



23_ (r) - 23_ (0) = In 
23+ (r) - 23+ (0) = [23_ (r) - 23_ (0)] - [23* (r) - 2?$ (0)] 



(t) = 7T / 
* Jo 



27r 2 T 



cot I i-^- + r ) 7rT + cot 



. 27ra; 



r 7rT 



(4.14a) 

(4.14b) 

(4.14c) 
(4.14d) 
(4.14e) 



where C « 0.577 is the Euler constant. 

To conclude this Subsection, let us proof the assump- 
tion of Sec. H about the Fermi liquid behavior of the 
system at low energies. In order to do this, we will cal- 
culate the fermionic Green functions, (^(t)^^)) and 

(V'p(t)V 'l(O) ), using the definitions flO]), (p~9|), and the 
results (4.14). Averaging over the bosonic fields similar 
to the well-known calculation of the Debye- Waller factor, 
yields: 



(4 (t)Vl (0)) = -L exp { ~ [23_ (r) - 23_ (0)] } X 
exp [ ~ [23+ (r) - 13+ (0)] j exp j ~ [23* (r) - 23* (0)] ] X 



™v{\[T>*+(t)+V*+(-t)] 



(4.15a) 



(^(r)^(0)) 



-^exp --[23+(r )+ 23+(Q ) ] x 
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cxp ji [V-(t) - V-(0)}\ exp |i [2?«(r) - Z>*(0)]} x 
e X p|i[^ + (r)-P $+ (-r)]|. (4.15b) 



Substituting Eqs. ( [4.14[ ) into Eq. ( f4.15| ), we obtain 



K(t)Vl(0)> 
(^(t)^l(O)) 



1 TFT 

27rw p sin 7rTr ' 



K(r): 



TXT 



sin 7rTr 



2ttvf 



exp 



■A-(r); 



(4.16a) 
(4.16b) 



efee coth (x+i7rrT) 



The Green function (4.16a) is not affected by interac- 
tions at all. The reason is, that it is taken at coin- 
ciding arguments (x\ — x 2 = 0), e.g., outside the in- 
teraction region. Because (^[(r)^i(0)) describes prop- 
agation of a chiral particle, the information about in- 
teraction is never carried back to th e observation point 
x = 0. The Green function ( 4.16bj ) acquires the free- 



fermion form at t > E c , which corresponds to the en 
ergies below the charging energy (note that T -C Eq)- 
In this energy range, (iP r (,t)iPl(0)) corresponds to a free 
fermion completely reflect ed fro m the dot. The phase 
factor exp(i2irAf) in Eq. ( 4.16b ) represents the scatter- 
ing phase 7r7V", which agrees with the Friedel sum rule 
(2.1). Thus, our intuitive picture of Sec. || is proven by 
explicit calculation of the fermionic propagators. 



B. Electrons with spin 

Similarly to the spinless case, we start here with the 
bosonization of electron operators: 

VR,a( x ) = CXp [ % 



'2-kX 



V2 



r] a I . y p L {x) + a<pj(x) 

lpL,a( X ) = ^^ CXp [ -l^ 



'2n\ 



V2 



(4.17) 



where index a — ±1 denotes the spin projections, and 
the Majorana fermions i]±i satisfy the anticommutation 
relations {77+1,77-1} = 0. Boson fields (p p L R and (p a L R 
corresponding to the charge and spin degrees of freedom 
respectively, satisfy the following commutation relations: 



-in sgn(x - y)<%; 



1 r{x), <pR(y) = in sgn(x - y)6ij\ 
l fiR(x),>fi J L (y) = -inSij] i,j = P, o. 



(4.18a) 
(4.18b) 
(4.18c) 



Like in the case of spinless fermions, it is convenient to in- 
troduce even and odd modes (p p ± for the charge and spin 



secto rs, and two x-independent fields, <I )P ' <T analogous to 
Eq. Q 



¥>l{x) = 7= ; i = P-,° 



V2 
V2 



(4.19) 



The commutation relations for the new fields wit hin the 
charge and spin sectors coincide with Eqs. (£h^); fields 
of different sectors commute wi th ea ch other. In terms 
of the new fields, Hamiltonian ( 3.15| ) acquires the form 
independent on dimensionless gate voltage J\f: 



Hn = 



<'F 



l=p,(T 7=± \ 



EE 



1 



(4.20) 



Backscattering Hamiltonian (3.17) takes the form 
2\r\v F 



Hbs = 



(0 P + (O) -izM) cos <pl(0), (4.21) 



and the effective action (3.16) can be rewritten as 
1 f p 

S = — dT 1 dT 2 L(r 1 -t 2 ) V 77 a (Ti)77 a (T 2 ) x 
nX Jo <£±i 

e 4 Q (^( Tl )-^(r 2 )) e 4a(# t '(r 1 )-4-(r 2 )) x 
e f(^(r 1 )-^(r 2 )) e f(i."(r 1 )-i."(r 2 )) x 



cos 



r ^(r 1 ) + ^(r 2 ) ^(n) + ^(r 2 ) 
h a 



sin 



-nJV - 
. (4.22) 



Similarly to Eq. ( 4.13 ), we introduce the relevant 
bosonic correlation functions 



V[_ (r) = (T T 0_(r)V_(O) i ); 
V\ (r) = (T T ^ + (r)^(0)); 
VI (r) = {T t ¥(t)&(0)); 
V% + (r)=(T r ^(r)^ + (0)), 



(4.23a) 
(4.23b) 
(4.23c) 
(4.23d) 



where index i = p,a labels charge and spin sectors re- 
spectively, bosonic fields are taken at x = 0, and av- 
eraging is performed over the Hamiltonian H$ given by 
Eq- (p0|). 

The calculation of these propagators can be performed 
immediately by n oticing that the spin sector of the 
Hamiltonian ( 4.20 ) corresponds t o the free bosons, and 
the charge sector differs from Eq. (4.10) only by replace- 
ment Ec — * 2Ec- Thus, we obtain 
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P£ to - Z?£ (0) = — / dare"* In 



£> p (r) - £> p (0) = (r) - (0) = In 



+ t) ttT sin (-if; +r) ttT 



A 



sinh' 
ttT 



( 7T 2 TX 



vf I sin7rTr 



V p (0) = In 



/ TTVp 



V p + (r) - V p + (0) = [X£ (r) - V p _ 
i r°° 7T 2 T 



[«M-r>»(o)] 



cot I i— + r J 7rT + cot I -i— + r ) ;, T 



T>1 = 0: 



> ( T ) = ysgnr, 



(4.24a) 

(4.24b) 

(4.24c) 
(4.24d) 
(4.24e) 



where C « 0.577 is the Euler constant. 

As we will see below, see also Ref. [14], the main con- 
tribution to the observable quantities is associated with 
the time scale r > \jEc- At this time scale the ef- 
fective theory can be further simplified. The mo de 



is "pinned" due to the charging energy, see Eq. (4.20) 
Therefore, the amplitude of quantum fluctuations of this 
mode is finite, see Eq. (4.24c), and the correlation func- 
tion (r) decr eases rap idly a t r > , as it fol- 
lows from Eqs. (|4.24d|) , (|4.24b| ) and fl4.24af ). The de- 
crease of correlations means that the average of a prod- 
uct, ( e ^+( Tl ) . . . e ^+( r «)^ can be replaced by the product 

of averages, (e l< ^+*- Tl ' ) ) . . . (e l( ^+' r "' ) ), if the intervals be- 
tween the times n, . . . , t„ exceed 1/ Ec- In other words, 



the operator functions of (p p + in Eq. ( 4.22 ) can be substi- 
tuted with c-numbers, according the rule: 



a " 4 



iZ>'(0) 



(4.25) 



On th e other hand, it follows from Eqs. (4.24a) and 
Q4.24b|) , that at r > l/E c 



(^(^M+^M-^^-d^O)) 2 ): 



2 In 



ttA / ttT V 
vfEcg c \smwTT J 



which means that such correlation function will be pre- 
served if we introduce another free bosonic field <p p (x), 
with commutation relation [(p p (x)] <p P (y)} — —iirsgii(x — 
y), and substitute 



backscattering term acquires the forml 







2\r\ ( E c e c v F 



1/2 



cos ttJV cos 0^(0), (4.28) 



and the effective action is given by 
1 [ p 

S = dT 1 d,T 2 L (n - r 2 ) V] Va(n)fi a (T2) x 

ir\ Jo a=±1 

e ia(^(r 1 )-^(r 2 )) e ia(* < '(Ti)-4 CT (r 2 )) e ^(Vp(n)-^(T 3 )) x 

^(ri)+^(r 2 ) 



a- 



-irAf 



tt ^(ri)-^(r 2 ) 
4 +a 2 



(4.29) 



Correlation functions of the free bosonic fields are given 
by 



V p {t)-V p {0)=VI (t)-V° ± (0)=ln 



A 



7TT 



Vf I sin7rTr| / 
(4.30) 



where cut-off A is of the order of vf/Eq because the 
charging energy Ec is the largest energy scal e whi ch can 
be considered with the help of Hamiltonian (4.27). It is 
easy to check also by an explicit calculation, that at time 
differences larger than Eq 1 correlation functions of the 
electron opera tors e valua ted with the help of the Hamil- 
tonians (4.20) and (4.27) respectively coincide. 



<pL(x = 0) + <l p -> V2<p p (x = 0). 



(4.26) 



Aft er substitutions ( |4.25| ) and ( [4.26] ), Hamiltonian 
(4.20) becomes just a Hamiltonian of three free bosonic 
fields 



dx 



d<P P 
dx 



( d(f+ 



\ dx 



d<p°L 
dx 



(4.27) 



V. THERMODYNAMICS OF THE "OPEN" DOT 



Coulomb blockade can be investigated experimentallyEZl 
by measu ring the differential capacitance of a dot, see 
Eq. ( PI) . In the regime of a developed blockade (weak 
tunneling between the dot and the electron reservoir), 
Cdiff{N) exhibits sharp peaks at half-integer values 
of Af. In the opposite limit of no backscattering, the 
differential capacitance is an A/"-independent constant, 
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Cdiffi-Af) = C. It was shown in Ref. [l4| that weak re- 
flection from a scatterer in the channel leads to the 
capacitance oscillations with a phase depending on the 
exact position of the scatterer. In this section we demon- 
strate that even at r — 0, the differential capacitance still 
depends on JV due to the electron backscatterings from 
inside the dot, which are described by the action S. The 
randomness of the backscattering events results in the 
randomness of the phase of the capacitance oscillations. 
We will relate the statistics of Cdiff{N) with the one of 
kernel L(t). 

The starti ng po int for the calculation of the capaci- 
tance is Eq. ( 3. IS ) for the thermodynamic potential. In 
principle, Eq. ( 3.18 ) enables one to consider the backscat- 
tering off a barrier in the channel (Hb s ^ 0), as well as 
off the dot (5^0). In the limit of weak backscattering, 
the perturbation theory in Hb s and S can be used to cal- 
culate C^ff{N). The case of S — was considered by 
MatveevEJ. He has shown that for spinless fermions, a 
non-vanishing result appears in the first-order perturba- 
tion theory, whereas for the spin- 1/2 electrons this order 
gives zero result. Similarly, it is sufficient to account for 
the scatterings from inside the cavity in the first order of 
the perturbation theory for spinless fermions, but in the 
case of spin-1/2 electrons we have to expand the thermo- 
dynamic potential ( 3.18| ) up to the second order in S, if 
Hbs = 0. 

Backscattering in the channel leads to a finite modu- 
lation of the average differential capacitance. The mod- 
ulati on am plitude can be estimatedtil by expansion of 
Eq. (3.18) to the second order in Hb s in the spin-1/2 
case. Electron scattering from inside the dot leads to the 
capacitance fluctuations superimposed on this modula- 
tion. The two contributions to the capacitance are not 
additive: the non-zero result appears only in the second- 
order in perturbations to the Hamiltonian Hq, when one 
expands Eq. ( 3.18| ) . In the domain of a relatively strong 
backscattering in the channel, |r| 2 > A/Ec, the leading 
term in fluctuations is proportional to the product of Hb s 
and S. We address the capacitance fluctuations at finite 
Irl in the end of this section. 



A. Spinless fermions 



1. Reflectionless contact 



The first-order expansion of Eq. ( 3.1 q ) in S yields: 

5n = T(T T S). (5.1) 



We substitute Eq. ( [4.12; ) into Eq. (|]|), retain on ly W - 
dependent terms, and obtain with the help of Eqs. (4.14) 



5n 



L_ [ dTL(j ) [e a ^K( T ) + c.c] , (5.2) 



where function K(t) is defined in Eq. ( 4.16b| ). To per- 
form the integration over r, we use the Lehmann repre- 
sentation (3.20) of the kernel L: 



sn = 



dt 



dr 



2-kvf J_ 00 2tt 
ttT 

sinh[7rT(i + ii 



(L R {t)-L A (t)) x 
^[e^K(r)+c.c.]. (5.3) 



The integration over r here can be now performed with 
the help of analy tic pr operties of function K(t). As it 
follows from Eq. ( 4.16b| ), the function K(t) is analytic in 
the lower semiplane Imr < 0. To calculate the integral 
of the first term in the brackets, we deform the contour 
of integration over r as shown in Fig. |[ 



P 



i t i) 



i t + p 



FIG. 3. The integration conto ur u sed in the evaluation of 
8Q, in the spinless case, see Eq. j5.3| ). Branch cuts of K(t) 
are shown by thick lines. 

Because of the periodicity of the integrand, the inte- 
grals over the parts of the contour running parallel to the 
imaginary axis, cancel out. As the result, only the pole 
contribution at r = it remains at t < 0. At t > the 
pole contribution disappears. Second term in Eq. ( |5.3| ) 
is integrated by using K{t)* = —K(—t) for th e real r. 
With the help of the explicit expression ( 4.16b ) for the 
function K(t), we find: 



i 



dt 



exp 

-.R 



2-kwf jo 

oo 

dxe~ 

o 



sinh(TrTi) 
„2tt 2 T 



(5.4) 



L H {t)e 



, / 2ttx 
coth I — — 

Eq \ E c 

L A {~ty 2 ^} . 



t ttT 



Finally, at low temperatures T <C Ec , the A/"-dependent 
correction takes the form: 



i 



dt 



2nivp jq 
[L R {t)e 2m * 



T exp 



L A (~t)e 



dx — 

o x ■ 

-27ri/Vl 



(5.5) 



Equation ( |5.5| ) relates the Coulomb blockade oscillation 
to the exact free-electron Green function in the dot. 
The variation of Sfl with the gate voltage is harmonic, 
however, its phase and amplitude are random quanti- 
ties. To reveal this oscillatory dependence in the aver- 
age quantities, one has to find the correlation function 
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5tt(J\f i)8&l(J\f2). At low temper ature s, Eqs. (5.5) and 
( 3.25 ) lead directly to the result ( 2.24 ) with the dimen- 
sionless function A# given by: 



(5.6) 



A(x) 



dy | 
-jexp 

V 



-xy + 2e y Ei(-y)] 



where Ei(x) = J^ oo e t dt/t is the exponential integral 



functi 



The correlation function of the differential 



capacitances (|3.19|) for different values of the gate volt- 
age and magnetic field is given by: 



8C dlff {l)5C dlff {2) _ 2 A 
C 2 E c 



A 



H 2 _ 



A 



Hi 

m 



cos 27m, 



(5.7) 



where we use the short hand notations i = Afi,Hi, 
n = Mi — A/2, and H± = Hi ± H 2 . Correlation mag- 
netic field H c is controlled by the charging energy and it 
is given by Eq. fl2.23| ). 

The variance of the capacitance fluctuations at H = 
is two times larger than in the unitary limit (H H c ). 
The dimensionless crossover function A(x) is plotted in 
Fig. |. 
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FIG. 4. Function A(x) determining correlation of differ- 
ential capacitances at different values of magnetic fields, see 
Eqs. Q and E q- (0)- 

Correlation between the capacitances at different mag- 
netic fields is suppressed if H- exceeds H c . In our ap- 
proximation, correlation "length" n c in the dimensionless 
gate voltage is infinite. To find n c , one should take into 
account that varying of the gate voltage affects the chem- 
ical potential of electrons in the dot by the level spacing, 
each time when J\f changes by one. If the chemical po- 
tential is shifted by Ec, a completely new set of levels 



determines the kernel L, thus suppressing the correla- 
tions. It results in the estimate: n c ~ Eq/A. Another 
effect which leads to the decrease of the correlation func- 
tion at large n, is the variation of the dot shape with the 
gate voltage. 



2. Finite reflection in the contact 

At fi nite b ackscattering in the contact, r ^ 0, Hamil- 
tonian (4.11) should be taken into account: 



Sft(N) = (Hi 



bs/- 



Calculatin g thi s average with the help of Eq. (4.14c), and 
using Eq. ( |3.19| ), we arrive to the result for the oscillating 
contribution of the averaged capacitance, 



SC(M) 

c 



2e l 



(5.8) 



This result was first obtained in Ref. F 



Average of the action Eq. (5T) simply adds to this 
result. Thus, we conclude that finite reflection in this 
order of perturbation theory does not affect t he m eso- 
scopic fluctuations of the capacitance, see Eq. (5.7). We 
will see below that this is qualitatively different from the 
electrons with spin where finite reflection in the contact 
leads to the increase of the mesoscopic fluctuations of the 
differential capacitance. 



B. Electrons with spin 

1. Reflectionless contact 

In the case of spin-1/2 fermions, the charging energy 
pins only one out of four m odes. Fluctuations spin mode 
are not suppr essed , see Eq. ( |4.27| ). These fluctuations av- 
erage action ( |1.29| ) to zero. To obtain a finite result for 
the A/"-dependent part of the thermodynamic potential, 
we have to expand Q up to the second order in 5*: 



5n = ——T(T T S 



(5.9) 



Upon the substituti on o f Eq. ([4.29[) into Eq. ((5.9|), we 
use the expressions (4.30) to perform the averaging over 
Hq . This cumbersome albeit straightforward calculation 
yields for the A/"-dependent part of the thermodynamic 
potential: 



50= 



{2-kvf) 



1 f 13 

Vf) 2 Jo ±ii 



ttT 



[sin7rT(ri + 170) sin7rT(r2 + ijO)] 
ttT 



1/2 



(5.10) 



[sin7rT(n + r 3 + ^0) sin7rT(T 2 — r 3 + «70)] 



1/2- 
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As in the Section V A , it i s convenient to use the 
Lehmann representation ( |3.20 ) for the kernel L: 



dti dt 2 



(2ttv f ) 2 J_ oa 2tt 2tt 
(L«(ti)-L A (ti)) (L R (t 2 )-L A (t 2 )) x 



(5.11) 



drx dT 2 dr 3 . , 
o smn7TJ (ti 



TFT 



7=±1 



iri) sinh7rT(i2 
ttT 



IT2) 



[sin7rT(ri + vyO) sin7rT(r2 + ijO)] 1 ^ 2 
ttT 



[sin7rT(ri + T3 + 270) sin7rT(r2 — T3 + £70)] 



1/2- 



The evaluation of the integrals over n and T2 here is sim- 
ilar to the procedure employed for the evaluation of the 
integral over r in Section V A . Upon the integration, we 
find: 



1 



(2ttv f ) 2 



dti dti 
2tt 2vr 



(5.12) 



(L R {ti)L R (t 2 )e- l2 * M 
ttT 

[sinh(7rTti)sinh(7rrt 2 

ft 

dr 3 

[sinh7rT(£i 



+ L A {t 1 )L A {t 2 )e^ N ) 



,1/2 



ttT 



ir 3 ) sinh7rT(t 2 + IT3)] 



1/2 ' 



Integral over t 3 can be easily evaluated with the help of 
the formula 



[sinh(x — itp) sinh(y + i<f>)] 



1/2 



\x+y\ 



where, K(fc) is the complete elliptic integral of the first 
kindba. We find for 



1 



(ttvf) 2 



dt\dti ttT 



[sinh(7rT£i) sinh(7rTi2, 



11/2 



(5.13) 



[L R (ti)L R (t 2 ) e - 12 ^ + c.c] e-^+'^K 



-7rr(ti+t 2 ) 



In the derivation of Eq. ( 5.12 ), we hav e utilized the re- 
duced version of action, see Eq. ( 4.29] ), which is valid 
only on a relatively long time scale larger than 1/Ec- 
Now, we average the product of two t hermodynamic po- 
tentials with the help of Eq. ( p\25| ). As we will see 
shortly, the result of the averaging is logarithmically di- 
vergent at large t. The divergence should be cut-off at 
ti,2 ^ l/^c- Without violating the logarithmical accu- 
racy, we can use also asymptotic expansion of the elliptic 
integral K(fc) = — In yl — k 2 . We obtain 

2A 2 

5fi(l)5fi(2) = t^TT COS 27171 X 



(2tt) 6 
dtidt 2 ir 2 T 2 
i/e c sinh7rT<i sinh7rTt 2 



In 



g-fl/TD +e -tx/T C 



T 2 (h+t 2 ) 2 



-til TO 



(5.14) 



Here, as in Eq. (5.7), we use the short hand notation 
i = Hi,J\fi, and n — Afi—J\f 2 . The diffuson and Cooperon 
decay times td and t q are related to the magneti c field 
values Hi, H 2 by Eq. (|3.24|) . Now we can use Eq. ( |3.19| ) 
to find the correlation function of the differential capaci- 
tances. In the leading logarithmic approximation we ob- 
tain 



SC cU ff(l)6C dt ff{2) 



E 

7=± 



C 2 



1-4 



8 A 2 4 (E c 
— k — t- hi -=- c°s 27m x 
3tt 2 E 2 c \T 



In max 1 



HI 



\ 



In(^) 



(5.15) 



where the fields H± = Hi ± H 2 are assumed to be sma ll 
compared to the correlation field H c given by Eq. (2.23). 
The new temperature dependent correlation field appear- 
ing in Eq. ( 5.1 5| ) is given by 



£0 / T 
A V 2nE T 



(5.16) 



with $0j -A, and Et being the flux quantum, the geo- 
metrical area of the dot, and the Thouless energy respec- 
tively. Unlike the scale H c , the characteristic field Hf is 
independent on the charging energy. This smaller field 
scale appears due to the existence of the "free" excitation 
mode <p+ which is not pinned by the effect of charging. 

With the increase of the number of channels in the dot- 
lead junction, the number of free modes also increases. 
The role of charging (which still pins only one mode), 
and therefore of the field H c , in the correlation functions 
should vanish gradually. The dependence of the correla- 
tion functions on the magnetic field at H < H c becomes 
power law rather than logarithmic; however, this power 
law is still non-trivial, and approaches Fermi-liquid re- 
sults only in the limit of an infinite number of channels. 

If H± > H c , [the charging correlation field H c is given 
by Eq. ( |2.23j )1, the correlation function of the fluctua- 
tion starts decreasing much faster, oc l/-ff±- Therefore, 
in order to get the representative statistics of the ca- 
pacitance fluctuations, averaging should be performed in 
the interval of magnetic fields larger than the magnetic 
field determined by the charging energy (2.23). Finally, 
in the limit of a strong field, Hi = H 2 > H c (unitary 

limit), the variance SC 2 ^ of the differential capacitance 
becomes four times smaller than at H = 0. 

Our results for the correlation functions diverge loga- 
rithmically at T — ► 0. At lower temperatures the pinning 
of the spin mode described by action-S* should be taken 
into account. A variational estimates shows that T at 
low temperatures should be replaced by Aln(£'c/A) in 
the above results for Hf and for the correlation function. 
We will elaborate on this point more in the end of the 
following subsection. 
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2. Finite reflection in the contact 

The main effect of the backscattering in the channel is 
that the Coulomb-blockade appears already in the aver- 
aged capacitances. Taking into account backscattering 
Hamiltonian (4.28) in the seco nd order perturbation the- 
ory, we obtain from Eq. ( [3.19] ): 



5C(J\f) = — |r| 2 cos27rAlnf — 
7T \ T 



(5.17) 



Due to the finite level spacing A, this result acquires 
mesoscopic fluctuations. As we already mentioned in the 
introduction to this section, at r ^ the leading term in 
fluctuations of the thermodynamic potential is first-order 
in both H{, s and S, 



6^ = -T(T T / dr Hts (t)S) . 
Jo 



(5.18) 



To calculate the average over the unperturbed state, we 
use the bosonized representation of Hb s and S given by 
Eqs. ( 4.28 ) and (4.29) respectively, and then the expres- 
sions (4.30) for the correlation functions of the boson 
fields. Then, similar to the Sect ion VA, we switch to 
the Lehmann representation ( [3.2C ) for the kernel L(t) to 
obtain: 



3(7rA0 



dt 



,2tt 

dridr 2 (ttT) 2 



(L R (t)-L A (t))x 



o sinh[7rTi + iirTin - r 2 )] [sin(7rTri) sin(7rTT 2 )] 1/2 
(ttT) 1 / 2 



E — 

±1 [sm wT(t 2 - n + ryO)] 



1/2' 



(5.19) 



Here only the A/"-dependent part of the thermodynamic 
pote ntial is taken into account. Integral over t± in 
Eq. (5.19) is determined by the contribution of the pole 
at n = t 2 — it, which can be easily calculated: 



dt {ttT) 1 ' 2 



n'y/TVVF ' J ^ [_j sinh(TrTt)] 1/2 

L A {t y^-l/4) + L RQ e -i«Vf-l/A) x 







dr-2 



TXT 



[sin(7rTr 2 ) sin(7rTr 2 + it)] 



1/2- 



(5.20) 



Integration over r 2 now is completely similar to the one 
we performed over the variable T3 in Eq. (5.12), and we 
find for 5 ft: 



■ cos(7rA/") 



(ttT) 1 / 2 \n{l/T 2 t 2 )dt 



[L A (t)e mM + L R (t)e 



i sinh(TrTt)] 



1/2 



(5.21) 



From Eq. (|^), with the help of Eq. Q3.25D and ( j3.19[ ), 
we find the correlation function of mesoscopic fluctuation 
of the capacitances: 



8C m (l)SC diff (2) _ 32e L |r| 2 A ^ ( E c 



C 2 



37T 2 Er 



T 



E 

7=± 



In max 1 



h: 



In 



V 



' Eg 

, T 



2X \ 



cos 27m x 



(5.22) 



where i = Hi,Mi, n = Af\ — A/2, an d H± = i?i±i? 2 . Cor- 
relation field H< is defined in Eq. (5. IE). The amplitude 



of fluctuations at a partial transmission of the channel 
is parametrically larger than at r = 0, cf Eq. (5.15). 
Furthermore, in the unitary limit the variance of the dif- 
ferential capacitance is suppressed only by a half of its 
zero-field value. This similarity with the case of spin- 
less fermions is due to the backscattering in the channel, 
which lea ds to pinning of the spin mode. 

Result ( 5.22| ) is valid at [Relatively high temperatures. 
As was shown by MatveevtJ, the divergences should be 
cut at energy e* ~ M 2 i?c cos 2 7rA/" which corresponds to 
the pinning energy of the spin mode. The higher-order 
corrections in backscatteringEJ shows that at T < | r | 2 
the logarithmic growth of fluctuations saturates. Si- 
multaneously, the correlation functions start to depend 
not only on the difference A/i — A/2, but also on each 
of these arguments separately. This weaker logarithmic 
dependence is beyond the scope of this paper. For an 
estimate of the differential capacitance variance at low 
temperature, one may replace \n(Ec/T) by ln(l/|r| 2 ) in 



Eq. J5.22D . 

Finally, we elaborate on the estimate of the character- 
istic energy scale which controls the low-temperatur e cut - 
off for the rcflc ctionless c onta ct. Compar ing E q. ( 5.17 ) 
with Eq. (|5.2lD, or Eq. (|5.21[) with Eq. fl5T5|), we ob- 



serve that the reflectionless case results can be obtained 
from formulas with finite reflection coefficient by putting 
|r| 2 ^ (A/Ec) ln(i?cA). It would correspond to the en- 
ergy of pinning of the spin mode e* 3^-A In {Ec / A) which 
agrees with our variational estimateE2l. 



VI. TUNNELING CONDUCTANCE OF THE 
"OPEN" DOT 

In the previous Section, we considered in detail the 
thermodynamics of the dot with one almost open chan- 
nel, and studied mesoscopic effects related to the dis- 
creteness of the charge. However, the majority of experi- 
mental work deals not with thermodynamics, but rather 
with transport through a dot. Coulomb blockade shows 
up as an oscillatory gate voltage dependence of the con- 
ductance of the dot connected with two leads. 

The case of small transparency of th^chaipels connect- 
ing the dot with leads is well studieaH'oliSO. The con- 
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ductance in the valleys can be represented as the sum of 
two physically different contributions - elastic and inelas- 
tic cotunnelingty. During the elastic cotunneling process, 
an electron enters the dot, spends there a time ~ h/Ec 
and then leaves the dot. Without interaction, this elec- 
tron would be able to spend a time of the order of U/A. 
As the result, the conductance in the Coulomb blockade 
valleys is suppressed by a factor of A/Eq- During the in- 
elastic cotunneling processES, an electron enters the dot, 
spends there time ~ h/Ec and then another electron 
leaves the dot. In this case, the final state contains an 
extra, in comparison with the initial state, two-particle 
excitation. It means that the phase volume of the final 
state is small as T 2 ; therefore, the elastic contribution 
dominates at T < (E C A) 1/2 . 




G 



R 



FIG. 5. Schematic view of the asymmetric two-terminal 
setup. The left point contact has one channel almost open 
and the conductance of the right point contact Gr is much 
smaller than e 2 /2Tvh. One of the electron trajectories con- 
tributing to the elastic cotunneling is also shown. 

Conductance of the dot connected to each reservoir by 
almost perfect channel was studied, in the limit A/ Eq — » 
0, by Furusaki and MatveevE3. They concluded that even 
small reflection in any of the channels leads to a dramatic 
suppression of the conductance of the system. At any 
gate voltage, except the discrete points of charge degen- 
eracy, they obtained a T 2 law which closely resembles 
the behavior of inelastic cotunneling in the weak tunnel- 
ing regime. Pursuing the analogy with the weak tun- 
neling regime further, it is natural to expect that there 
should be another non-vanishing at T — > contribution 
from a counterpart of the elastic cotunneling mechanism. 
Studying of this contribution is the goal of this sectioncil. 

Similar to Ref. two limiting cases may be distin- 
guished. In the first case both channels are either open or 
have the same reflection coefficient. This case is techni- 
cally difficult to consider. Instead, we concentrate on the 
properties of the strongly asymmetric setup; one point 
contact (left in Fig. |J) has the transparency close to unity 
and the other contact (right in Fig. ||) has a very small 
conductance Gr, <C e 2 /irh. This case can be realized ex- 
perimentally by a corresponding adjustment of the volt- 



ages on the gates forming point contacts^. Moreover, it 
follows from the scaling argumentsEII that the asymmetry 
of contacts is a relevant perturbation. In the case A = 0, 
the strongly asymmetric limit corresponds to the fixed 
point of a system with an infinitcsimally small initial 
asymmetry. Therefore we can expect that at A/Ec <C 1, 
this limit adequately describes dots with a finite initial 
degree of asymmetry. 



A. General formalism 

For calculation of such tunneling conductance, we have 
to modify derivation of Sees. Ill A and IIIB in order to 



take into account the tunneling between the dot an d the 
second lead. In comparison with Hamiltonian (3.1), the 
total Hamiltonian of the system acquires two additional 
terms, 



H = Hp + He + H 



M 



H 7 



(6.1) 



Here Hp describes the electron motion in the dot and in 
the left lead and is given by (^^), interaction Hamilto- 
nian He is given by Eq. (1.1), and Hm is the Hamiltonian 
of free electrons in the right lead 



P 



(6.2) 



Tunneling Hamiltonian Ht describes the weak coupling 
between the right lead and the dot, 



H T = v t 4^(r t )'Y]a p + h.c. 



(6.3) 



where r t is the coordinate of the tunneling contact, ^(r) 
is the electron wavefunction in the lead, and Vt is the cou- 
pling constant which will be later related to the tunneling 
conductance of the contact Gr. 

Because Gr <C e 2 /(2ttK), we can consider the tunnel- 
ing current I as the function of applied voltage V in the 
second orde r of perturbation theory in tunneling Hamil- 
tonian ( |6.3|) . This gives us the standard resulto 

J(eV) = i [J («Q„ -> eV + iO) - J (iO n eV - i0)\ , 

(6.4) 

where Q n = 2irTn is the bosonic Matsubara frequency, 
and Matsubara current J is defined as 



j(»n„) = ev\v / dTe- in ™ T g M (T)n(T). 

Jo 



(6.5) 



Here v is the one-electron density of states per unit area 
and per one spin in the dot, Qm is the Green function of 
the electrons in the leads, 



g M = -J2(T T a Pl {r)al(0)) 

P1-.P2 



ttT 



sin7rTr 



(6.6) 
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with vu being the one-electron density of states per one 
spin in the lead, and function II (r) is given by 



n(r)=z,- 1 (T r ^(r;r t )V(0;r t )). 



(6.7) 



Averages in Eqs. @ and (|^) are performed with re- 
spect to the equilibrium distribution of the system with- 
out tunneling. We choo se to introduce v into Eq. (6.4) 
and into definition (6.6) to make function II(r) dimcn- 
sionless. 

In the absence of the interaction, Ec = 0, propagator 
n(r) is nothing but the Green function of non- interacting 
syste m; its ensemble average has the form analogous to 
Eq. ©, 



n(r) 



nT 



e c =o sin7rTr 



(6.8) 



Then, substitution of Eqs. fcj) and (3.8) into Eq. fl6.5|) 



and analytic continuation (6.4) gives the tunneling cur- 
rent / = sGrV (s is the spin degeneracy), where the 
tunneling conductance of the contact per one spin is 



Gr 



2ire 2 



-v t v M v- 



(6.9) 



With the help of Eq. ( |6.9| ) we can rewrite Eq. (|6_.5| 
terms of the bare conductance of the point contact 



J(ifl r , 



Gr 
2ire 



o 



, 7rTe- l °" T , , 

dr— n(r). 

sin itTt 



(6.10) 



As we will see below, function Il(r) can be analyti- 
cally continued from the real axis to the complex plane, 
so that the result is analytic in a strip < Re r < (3, and 
has branch cuts along Rer = 0, (3 lines. It allows one to 
deform the contour of integration as shown in Fig. [g, and 
to obtain 



J{iVL n ) 



GrT 
2e 

/ U(it ■ 



die""' 



0) 



e{-n n )d(t) - 6(n n )9{-t)} x 

U(it - 0) 



Vsinh[7rT(t - i0)] sinh[7rT(f + iO)] 



(6.11) 



Now the analytic continuation ( |6.4| ) can be performed, 
because the periodicity of the Matsubara Green functions 
was already taken into account. This gives 



I(eV) = i 



GrT 
2e 



dtc 



-ieVt 



IL(it + 0) 



IT(it-O) 



sinh[7rT(t - i0)\ sinh[7rT(t + i0)} 



.12) 



Next, we use the analyticity of II(r) in the strip < 
Rer < (3, and shift the integration vari able t — > t — i/3/2 
in the first term in brackets in Eq. (5.12), and t — > 
t + iP/2 in the second term. Bearing in mind that 
n(r) = -n(r + j9), we find 



/ = Tsinh 



eV 
2T 



Gr 



die 



-ieVt 



n ( u + 4 



cosh nTt 



.13) 



Linear conductance G is therefore given by 

n(it + f) 



G — Gj 



dt- 



a) 



2cosh?rTi 
b) 



(6.14) 







FIG. 6. The integration cont our used in the evaluation of 
the conductance, see Eq. (5.10) for (a) f2 n < 0, and (b) for 
fl n > 0. Branch cuts of the analytic continuation of II (r) are 
shown by thick lines. 

Let us turn now to the actual calculation of the func- 
tion II(t). It was shown in Ref. E7] that the interaction 



drastically affects the form of the function (6/7), however, 
some contributions were not taken into account. Our 
purpose is to construct an effective action theory, similar 
to that of Sec. Ill, for calculation of LT(r). Once again, 
we wish to get rid of the fcrmionic degrees of freedom of 
the dot. Similar to Eq. (3.4), it is convenient to rewrite 



charge operator in terms of the variables of the channel. 
However, here we have to keep in mind the fact that the 
tunneling events described by operators ip' (r t ) and ip (r t ) 
change the charge in the system by +e and — e. It can 
be taken into accounted by introducing three additional 
operators: Hermitian operator n, and unitary operators 
F, with the following commutation relations: 



n,F 



F 



(6.15) 



We can definitely choose Hilbert subspace in a way such 
that operator n has integer eigenvalues. Finally, these 
operators commute with all the fermionic degrees of free- 
dom. Then, we can change the definition of the charge 
operator [cf. Eq. (fO)] to 



Q 



channel 



drip^ip + n, 



(6.16) 



and rewrite Eq. (6.7) as 



n(r) = v-\T t F{t)^ (r; r t ) F(0)V (0; r t )>. (6.17) 



It is easy to see from Eqs. Q6.16 ) and ( |6.15 ) that opera- 
tors P\F in Eq. ( |6.17 ) change the charge by +e and — e 
respectively, in accordance with the initial definition of 
charge. 

After this manipulation, the Hamiltonian of the sys- 
tem and correlation function ( |6.17 ) become quadratic in 
the fermionic operators of the dot, so that part of the sys- 
tem ca n be integrated out. We use the identity similar 
to Eq. (|3~6|): 
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; {e-^T r ^ 2 (r)V 2 (0)} 
e"^ 2 (T T ^ 2 (r)^ 2 (0)) 2 e-' 3 ^T T e^ -C^ 1 /><^(^)^^» = 



(6.18) 



e -/3n 2e -/3H! 



(It?, I <ItaT^ 



(r^ 2 (r)^ 12 (r3)> 2 (T T ^ 2 (0)^ 12 (r 3 ))) 2e */o Vl / l >^^^)^^)> 2 , 



n 



where all ip 2 ,tp2 are the fermionic operators of the dot, 
and the rest of the notation is the same as that in 



Eq. (|3_6|). 

The calculation of the product (ffi 2 (r 1 )ff 12 (r 2 )) 2 was 
performed in Sec. Ill A, see Eq. (3.7), and all the s teps 



leading to the derivation of the effective action ( |3.14[) can 
be repeated here. Calculation of the remaining operator 
products can be performed along the lines of Appendix [a|. 
This yields 

(T T V(ri;r i )^ 12 (r 2 )) 2 = -^{t 2 ,Q)R* ( Tl - r 2 ), (6.19) 
(T r ^(ri;r t )^ 12 (T 2 )>2 = -^(r 2 ,0)i2(r a - t x ), 



THl ip(x)e~ rHl are 



where, similar to Eq. (3.7), %1>{t\x) = e 
the one dimensional fermionic operators of the channel 
in the interaction representation, V>(r) = t). Ker- 
nel R describes the motion of an electron from the tunnel 
contact to the entrance of the single mode channel, and 
it is given by 



R(r) = f dy^y)d x G(T;r,v t ). 
2m J 



(6.20) 



Here, Q is the exact Matsubara Green function of the 
closed dot subjected to the zero boundary condition. The 
wave function <p(y) describes the transverse motion in the 
single-mode channel, and the coordinates x in the deriva- 
tive of the Green function Q is set to +0. 

Kernel R(j) is the random quantity with the zero av- 
erages. In the universal regime, products of retarded 
R R (t) and advanced R A (t) counterparts of R(t) entering 
into the Lehmann representation (3.20) have the follow- 
ing non- vanishing averages: 



(R% 1 (t 1 )R£ 9 (t 2 ))=Av F 8(ti+t2)0{t 1 )e t h, (6.21a) 



AujP$(ti-t 2 )0(ii)e t h, 



(6.21b) 



where the decay times t h ' D associa ted w ith applied mag- 
netic fields H± t 2 are given by Eq. All the higher 
momenta can be fo und by using the W ick t heoremE! 
Deriving Eq. (|6.21b|) we use Eqs. (|3.2l|) , ( |3.23| ), and the 
identity a R (t;n,r 2 ) = [G A {-t; r 2 , n)] . 



To complete the derivation of the effective theory, we 



use Eqs. (6.18) and (6.19), int roduc e left and right mov- 



ing fermions similarly to Sec. (Ill A), and thus obtain the 
effective action representation for n(r) from Eq. (3.17): 



II(t) =n i „(r)+n ei (r) 
n,„ = - 



G{ ~ T;Tt > rt) (T T e- S F(T)F(0)); 



v{T T e 



(6.22a) 
(6.22b) 



i r 

, g- / dTxdT 2 R{Ti-T)R*{-T 2 )x (6.22c) 

v{T T e b ) Jo 



(T T e- § F( T )F(0) [Mn) + i> R ( ri )] [^(r 2 ) + ^r(t 2 )]). 

Here the averaging is performed with respect to the 
Hamiltonian 



Ho = ivF dx 



{^Ml-^bP^r) (6.23) 



T 



+ ^ ( / dx: ij) ] L ^ L + ip R ipR :+M -n 



and action S is given by Eq. ( 3.16 ). The difference be- 
tween the Eq. ( 6.23 ) and Eq. ( [3. 15 ), is caused by the 
differ ent d efinitions of the charge operator in Eqs. (3.4) 
and ( 6.16 ). 

Two contrib ution s can be distinguished in the co rrela- 
Inelastic contribution ( 6.22b| ) was 
in the approximation correspond- 



tion function ( p. 22 ) 
considered in Rcf. ^ 
ing to S = 0, and with the Green function of the dot Q 
replaced by its averaged value Q. The obtained results 
vanish at low temperatures. The reason for the vanish- 
ing is that this term does not allow the introduced elec- 
tron to leave the dot; charge of the dot at the moment 
of tunneling suddenly changes by +e and all the other 
electrons have to redistribute themselves to accommo- 
date this charge. The logarithmical divergence of the 
imaginary time action corresponding to such evolution 
(orthogonality catastrophe) completely suppresses this 
contribution at T — > . Con versely, the second contri- 
bution, n e ; from Eq. ( |6.22cj ), contains the kernel R{t) 
which promotes an electron from the tunneling contact 
to the channel. Because the very same tunneling electron 
is introduced to and then removed from the dot, there is 
no need in the redistribution of other electrons, so no or- 
thogonality catastrophe occurs. As the result, the elastic 
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contribution survives at T — > 0, analogously to the elastic 
cotunneling contribution for the weak coupling regime. 

In what follows, we will be interested in the low tem- 
perature beh avior o f the system, so we will retain elastic 
contribution (6.22c) only. Similarly to the Sec. [Tv| results 
for electrons with spin and spinlcss electrons differ signif- 
icantly, and we will consider those two cases separately. 



exp 



1 7= 

y/2 4 



exp 



.$(ti) - $(r 2 ) 
1 7= 

0+(t 2 ) , 7T 



V2 



We will consider only the elastic contribution ( 6.26b| ) 
because it does not vanish at low temperatures. 



B. Spinless electrons 



We follow the lines of Sec. IV A in the bosonization of 



the chiral fermionic fields. In order to account for the 
appearanc e of t he operato r h in the Hamiltonian [com- 
pare Eqs. ( 3.1 5| ) wit h Eq. ( 6.23 )1, we change slightly the 
transformation (4.9): 



<Pl{%) — t= wv + 7rn; 



V2 
V2 



(6.24) 



where operator n commutes with the bosonic fields 
<p±,<f>. In order to preserve the commutation relation 

F j <Pl,r = 0, we change the operator as 



Ft »&e- iy/ **,F>->Fe< 



V2& 



(6.25) 



The fact that c omm utes with boso nic fields <pL,R is 



obvious from Eqs. ( 6.15 ) and Eq. ( 4.8d ) 



Su bstitut ion of Eqs. Q), ( |6.24| ) and ( |6.25| ) into 
Eqs. ( |6.22b| ) yields 



n in (r) = - 



g(-r;r t ,r t ) 
v{T T e~ s ) 



(T T e- 5 F(r)F(0)e^[*^-*( T )]) 



(6.26a) 



for the inelastic part of the cotunneling, see also Ref. |27 
For the elastic contribution, we find 



n ei (r): 



?Kn)^2)(-ir (Tl)+fi(r2) 



l — s - ( dndT 2 R{n - t)R*(-t 2 ) x 
Ve b ) Jo 

i+4>(ti)-4(t 3 )] x 
ip-{n) - <p-{t 2 )' 



(T r e-^(r)F(0)e^[ 2 * (0) - 2 * (r)+ * (Tl) -* (T2) ] x 



exp 



V2 



n 



cos 



y/2 4 



(6.26b) 



Averaging in Eqs. (6^26). is performed over the Hamil- 
tonian give n by Eq. ( 4.1C ) , fi nite backscattering is de- 
scribed by (4.11), and action (4.12) is modified as 



S=^- fdT 1 dr 2 L{r 1 ~T 2 )f 1 {T 1 )f 1 {T 2 )(~lf^ +il ^ 
2ttA Jo 



1. Reflectionless contact 



In the lowest in A/Ec approximation we can ne- 
glect action S in Eq. ( |6.26b| ) at all. Then averaging 
over boson ic fields c an be performed with the help of 
Eqs. (|OJ) and Eq. ( gig ), average of the product Ma- 
jorana fermions operator is given by Eq. ( |4.7| ), and the 
relevant correlation function of the operators h, F, F^ is 
given by 



(T T F(T)F(0)(-lf^+^) = sgn(r - n)sgn(T - r 2 ). 

(6.27) 



Equation ( |6.27 ) follows from Eq. ( 3.15 ) and from the fact 
that operators n, F commute with Hamiltonian (4.11) 
and thus do not have their own dynamics. We obtain 



n e; (r) 



2tt\K(t)\ 2 
nT 



rP 

/ dT X dT 2 R(Ti - t)R*(-t 2 ) x 
Jo 

K( T2 )K( n - t) 



sin7rT(ri — t 2 ) K(t\)K(t 2 — r) 



(6.28) 



K( n - t 2 ) 



K{-t 2 )K( Ti - t) 
K{ n )K(r-r 2 ) 



+ c.c. 



where function K(r) is defined by Eq. ( 4. 161 ), and C is 
the Euler constant. 



a) 



o x 



b) 



i L+ p (• 
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X T 2 p 

i t A + p 



1 



FIG. 7. The integration contour used in the evaluation of 
the conductance in the s pinle ss case, for the first (a) and the 
second term (b) in Eq. ( |6.29| ) . Branch cuts of the integrand 
are shown by thick lines. Contribution of the semi-pole at 
Ti = T2 in the first term in brackets in Eq. (3.29) is canceled 
by the term complex conjugated to it. 
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Be fore performing the analytic continuation, see 
Eq. ( 6.14 ), we hav e to transform integrals over imagi- 
nary times in Eq. ( |6.28 ) to the integrals over real time . 
In order to do so, we use Lehmann representation ( 3.2C ) 
for the kernel R(t): 



n ei (r) = 



2tt\K{t) 



dti 
1^ 



2tt 



(6.29) 



R A (t 1 )) (R R (t 2 )~R A (t 2 ))' 
n 2 T 2 



(IT l d,T 2 — — — — 

i) sinh|7rT(£i + in — ir)J sinh[7rT(t2 + iT 2 )\ 



7rT 



K(t 2 )K( Ti - t) 



sin7rT(ri — t 2 ) K(t\)K(t 2 — r) 



K(ti - t 2 ) 



K(-t 2 )K( Ti - t) 

K( Ti )K(t-t 2 ) 



+ c.c. 



Int egrat ion can be now performed in a manner similar to 
Sec.VA. Using the fact that function K(ji) is analytical 
within the lower complex semiplane Imri < 0, we deform 
the contour of integration as shown in Fig. ^. 

Because of the periodicity of the integrand, the inte- 
grals over the parts of the contour running parallel to 
the imaginary axis, cancel out. As the result, only the 
pole contribution at n = it± + r remains at t% < 0. At 
t\ > the pole contribution disappears. Analogously, 
the complex conjugated terms are contributed by pole 
Ti = it i + t at <i > 0. As the result, we obtain from 
Eq. ( p9D 



n ei (r) 



2w\K(t)\ 2 f dtxdh (uR 
2tt 



ttT 



ttT 



(R R (t 2 )~R A (t 2 )) / dr 2 - . 

( R A (t 1 )K(r 2 )K(it 1 ) R R (t 1 )K(-r 2 )K(-it 1 ) 
sin7rT(T - r 2 + if x ) + t)K(t 2 - r) K(-it x - t)K(t - t 2 ) 



(6.30) 



K(t -T 2 + itl 



| R A {—t 2 )K (iti) _ e _ i2nAr 

K(itl+T)K(T-T 2 ) 



K(T2-T-Ui) 



R R (t 1 )K(T 2 )K{-it 1 ) 
K{-it x - t)K{t 2 - t) 



where we wrote explicit expressions for all the terms. Integration over t 2 can be now easily performed by deformation 
of the integration contours shown in Fig. 0, and we obtain 



n e; (r) 

2it\K(t)\' 



WFE},e 2C 



2it\K(t)\ 2 
vvpE^e 20 

dt\dt 2 



f dt 1 dt 2 {R R (t 1 )[R R (t 2 )]* + R A (t 1 )[R A (t 2 )}*} 



TTT 



ttT 



sin ttT(t — it 2 + it\) 



Kir ih)[R W^HW*) 



T(it 2 - iti - t)] 

R A (t l )[R A (t 2 )}*K(it 2 )K(it 1 ) t RXjhftRXfoWKi-itdKj-ih) 

K(-it x - t)K(t - it 2 ) 
,R R (h)[R A (t 2 )]*K(it 2 )K(-it 1 ) 



(6.31) 



K(it!+T)K(T -it 2 ) 



K(it X +T)K(it 2 -T) 

~ l27TM K{it 2 - it t 



K(-it 1 - r)K(it 2 - t) 
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f (T t 
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FIG. 8. The integration c ontou rs used in the evaluation of 
the integral over T2 in Eq. (6.30) for (a) the first and third 
terms in brackets, and (b) for the second and fourth term. 
Only first and second term have poles at T2 = r + it\. 



One can see from Eq. ( |6 .31 ) that function H e i(r) is in- 
deed analytic in the stri p < Rer < (3, which justifies 
the steps leading to Eq. ( |6.14 ). 



Finally, we substitute Eq. fl6.3l| ) into Eq. ( |6T4). Fo r 
small temperatures T -C Ec, we have from Eq. ( 4.16bj ) 
K ((3/2 + it) = txT/ coshnTt. As the result, the first term 
in Eq. ( |6.31 ) produces a contribution cx T 2 and can be 
neglected. The remainder can be recast into the formula 



G = 



2itG r 
WFE"l,e 2C 



/ dtK(it) [R A (t)e l2 ^ + R R (~t)} 

J — OC 



(6.32) 



which gives non-averaged conductance of the dot. Here 
we used the fact that the characteristic scale of the in- 
tegration over ti, t 2 ~ 1/Ec is much smaller than (3. 
Equation (6.32) is reminiscent of the Landauer formula. 
However, the form-factor K(it) entering into this formula 
indicates that a large number of states in the dot partic- 
ipate in the transport, unlike the case of non- interacting 
electrons. 
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Now, we are prepared to study the statis tics of the 
conductance. Using the expli cit expression ( 4. 16b|) for 
function K and formula ( 6.211: ), we find the average con- 
ductance 



_ 2A 
G = G R —e 
bj C 



2C 



A(0), 



(6.33) 



where C w 0.577 is the Euler constant, and A(0) w 1.398 
is given by Eq. ( |5.6| ). This expression is analogous, to, the 
elastic cotunneling for the case of weak couplingi3'EJ. 

For the correlation function of the mesoscopic fluc- 
tuati ons o f the conductance, we find with the help of 
Eqs. fl6.2lD 

2 r 



SG{1)SG(2) 



G 2 



cos nn 



A 



Hi 



A 



Hi 



(6.34) 

where we use again the short hand notations i = Mi, Hi, 
n = Mi — M 2 , and H± = H\ ± H 2 . Correlation mag- 
netic field H c is defined in Eq. (2.23), and the dimen- 



sionless function A(x) is given by Eq. (5.6) and is plotted 
in Fig. |]. Once again, we see that even though the aver- 
aged conductance does not any longer oscillate with the 
gate voltage, the discreteness of charge manifests itself 
in the oscillatory behavior of the conductance correla- 
tion function. It is also noteworthy that the mesoscopic 



conductance fluctuations are of the onier of the average, 
similarly to the weak coupling regimetil. 



2. Finite reflection in the contact. 



So far, we have shown that the conductance in the tun- 
neling setup is non-vanishing at T — > 0, which is anal- 
ogous to the elastic cotunneling in the weak tunneling 
regime. However, the oscillatory dependence of the con- 
ductance showed up not in the average conductance but 
rather in the correlation function of mesoscopic fluctu- 
ations. On the other hand, as we saw in Sec. [v] finite 
backscattering leads to the oscillatory dependence in the 
averaged quantities. The purpose of this subsection is 
to study how the finite reflection in the contact affects 
the elastic cotunneling, and to demonstrate that it in- 
deed leads to the oscillatory dependence of the averaged 
conductance on the gate voltage. 

To treat the finite reflection in the contact, we have to 



expand the denominator and numerator in Eq. ( 6.261: ) 
up to the first order in the backscattering Hamiltonian 
( 4.1 1|) . Performing averaging o ver th e boso ni c fiel ds, we 
obtain with the help of Eqs. flP4| ), ( |4~S| ), ( |6~27| ), and 
( H.16bp : 



n ei (r) 



\K(i 



2-kvvfEcc c 
K{t 2 )K{ Ti - T 



/ dndT 2 dT 3 R{n-T)R*{-T 2 )- — - 
Jo sin it, 



X 



siu/iT(ti — T2) 

-^K(t 3 )K(t - r 3 ) sin7rT(r 2 - t 3 )K(t 2 - t 3 ) 

+ c.c. 



(6.35) 



[K( Tl )K(r 2 - t) \ K(-t 3 )K(t 3 - r) si n7 rT(ri - t 3 )K{t x 
e^K(-r 3 )K(r 3 - r) sin 7 rT(r 1 - t 3 )K{ Ti - r 3 ) 



K(t 3 )K(t - t 3 ) sin7rT(T 2 - t 3 )K{t 2 - r 3 ) 



- 2 cos 2ttM 



Here we retained only the terms which do not vanish after ensemble averaging. Then, we can use the Lehmann 
representation for the kernel R(t) and perform the integration over n 2 in the manner of the previous subsection. It 
yields 



n el (r)=-\r\ 



\K(t)Y 



2TrvvpEce c 



dtidt 2 



dn 



nT 



r R A (h)[R A (t 2 )]*K{it 2 )K(it 1 ) 
K(ih + r)K(it 2 - r) 

r 3 ) 



K{t 3 )K{t - T 3 ) 



sin7rT(r — it 2 + iti) 
sinTrT(it 2 — T 3 )K{it 2 



K{—t 3 )K{t 3 — t) sin7rT(r + it\ — t 3 )K{t + it\ — t 3 ) 



K(-t 3 )K(t 3 - t) sin7rT(r + it x - t 3 )K(t + it x - t 3 ) 



K{t 3 )K{t - T 3 ) 



smirT(it 2 — r 3 )K(it 2 — r 3 ) 



- 2 cos 2ttN > + c.c 



(6.36) 



(6.37) 



part of the ensemble averaged conductance 



We deform the integration contour over r 3 , a s shown 
in Fig. [|, and substitute the result into Eq. ( 6.14 ). Then, 
we average the p roduct of the Green functions R A [R A ]* 
using Eq. (6.21b). As the result, we find the oscillating 



G = ai G R 



Er 



i27r7V. 



Here ct\ is the numerical coefficient given by 



(6.38) 
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4e 



c 



Oil 



dxdy 



1 sin(-7re 



sinh [e~ v Ei{y) + e x+v Ei(-x -y)- e v Ei(-y)] w 1.458, 

with Ei(x) b eing the exponential-integral function@. 

Equation ( 6.38| ) confirms our expectation that the fi- 
nite backscattering leads to the oscillatory dependence of 
the averaged conductance on the g ate voltage. Although 
the amplitude of oscillations ( 6.38 ) is s mall compared to 
the average value of the conductance ( 3.33| ), it still ex- 
ceeds at low temperatures the contribution of the inelas- 
tic cotunnelingrH to the conductance oscillations. 



FIG. 9. The integration conto ur used in the evaluation of 
the integral over T3 in Eq. (3.36). Integral along the part of 



the contour parallel to the real axis cancels out the 2 cos 2nN 
term. 



C. Electrons with spin. 



We use formulas of Sec. IV B to bosonize the chiral 



fermionic fields. In order to account for the appearance 
of the ope rator h in the Hamiltonian [compare Eq s. (|3Tl5| ) 
with Eq. ( |6.23 )1, we change the transformation (|4.1S| ): 



(p\{x) + g_Qr) + - 6 ip ir(Af - h) 
71 ; 



ip R (x) = -j= , (6.39) 

where i = p,a, and operator n commutes with the 
bosonic fields (£ P t; ' cr , <3E> P,CT . To preserve the commutation 

relation F^,(p p £ a R — 0, we change the operators F',F 



similarly to Eq. ( 6.25 ) 



F+ h-» fU-^" ,F h+ Fe' 



<i><> 



(6.40) 



Subst itution of Eqs. (4.17), (6.39) and (6.40) into 



Eqs. (6.22b) results in the formulas similar to Eq. (3.26): 



n m (r) = - 2G [ T;rt l rt) (T Te ^F(r)F(0) e '[*(»)-*(-)]), 



v{T T e~ s ) 



(6.41a) 



n ei (r) 



g- [ dTidT 2 R(n-T)R*(-T 2 ) V (T T e- i l i (T)J , (0)^ a (ri)^a(72!)e i [*'W-*'W]] x 

r re- s ) Jo a=±1 



TTUX{T T 

e fQ(^(r 1 )-y-(r 2 )) e fa(4."(r 1 )-*-(r 2 )) e §( V 3 p _(r 1 )-^(r 2 )) e ±(^(r 1 )-*''(r 2 )) x 

^(n)-^(r 2 ) ^(ri)-^(r 2 ) tt 
2 + a g 2 ^ ^ ~ 



cos 



p + (ti)+0 p + (t 2 ) V+{n) + ^ + {t 2 ) tt 



+ 2 («( T i) +«(r 2 )) -ttAA 



(6.41b) 



for the inelastic and the elastic cotunneling respectively. 
However, for the calculation it is more convenient to pro- 
ceed directly to the low energy effective theory (4.28), 
because the main contribution to the conductance comes 
from the energy scale smaller than the charging energy 
Ec- First, we integrate out the symmetric charge mode 
(p p + . Then, we wish to use the substitution (4.26). The 
important difference brought into the problem by ac- 
counting for the second junction, is that Eqs. (6.41) con- 



tain the charged field $ p itself, and not only the combi- 
nation $ p — 4> p _{x = 0), as we had before. Fortunately, 
the corresponding change can be simply accounted for 
by the introduction of one more chiral field <fi, so that we 



have 



p _(x = 0) + l> p -> V2ip p (x = 0), 
<f> p - <p p _(x = 0) -> \/2<j>{x = 0) 



(6.42) 



instead of Eq. ( 4.26 ). The field <j> enters neither the ef- 
fective action, nor the backscattering Hamiltonian, so it 
can be immediately integrated out, and we find the low 
energy representation 



s i[$(x)-4(0)] ^ / e ilk-r)-k0)]/V2K e il^(r)-<p"(0)]/V2 



A V sin ttTt 



1/2 



(6.43) 
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The prefactor in the last formula can be found by re- fectiv e the ory a nd th e initial theory to coi ncide . Using 
quiring the averages calculated with the help of the ef- Eqs. ( 6.42 ) and ( 6.43 ), we obtain from Eq. ( 6.41 ) 



n m (r) 
n ei (r) 



2g(-r;r t ,r f ) n fv£ f ttT \ 1/2 s p ^ e i[^(r)-#'(o)]/^ 

1 / / \ 1/2 -a 

1 7T Vp / irl \ 



'{T T e- s ) E c e^ \ A 3 



— — - / dndT 2 R{n - t)r*(-t 2 ) V (r r e- s ^(r)^(o) 



(6.44a) 
(6.44b) 
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7T, g+(n) - y+( T 2) 7T ' 

4 + 2~^ + 2 ~ ^ 



where the averaging is performed with respect to the Hamiltonian (4.27). The effective action ( 4.29Q acquires the 
form 



§ = -1 [ d,T i dT 2 L ( n - T 2 ) V fjcintfcfo) X e 4«(^(n)-^ (r 3 )) e ia(^(n)-4<'(r 2 )) e ^(^(r I )-^(r 2 )) ^ ^ 

ttA Jo q=±1 



cos 



^(Tl) + ^(T 2 ) 7T ■ 

a + - (n(n) + n(T 2 )) - irJV 



cos 



- + a + - (n{Ti) - n(T 2 )) 



y§ 



and the backscattering Hamiltonian H\, s is given b; 



1/2 

//,,. = ^ ( (-if cosxAfcos^O), 



7T V TTA 



(6.46) 



cf. Eq. ( ]£2Sj ). Cut-off A in Eqs. - ( |g§ should be 

of the order of vp/Ec, because the charging energy Eq 
is the largest energy scal e whi ch can be considered with 
the help of Hamiltonian (4.27). 

In the f ollow ing subs ection s we will apply effective de- 
scription Q6.44D - Q6.46D and ( |6.14| ) to find the tunneling 
conductance of the dot in the asymmetric setup. 



1. Reflectionless contact. 



Let us consider first the elastic contribution ( 6.44b| ) 
(we will see below that inelastic contribution should be 
taken into account in order to obtain the correct temper- 
ature dependence). In the lowest-order approximation 
we n eglect the action S, and obtain with the help of 
Eqs. flO0[): 



n£?(T) = 



i 



ttT 



vvpEc& sin7rTr j 

R{n - t)R* (-t 2 ) 



dr\dr 2 x 

TTT 



(6.47) 



sin7rT(ri — t 2 ) 



/ sin ttTti sin ttT(t 2 — r + Z70) \ 
V sin7rTr2 sin7rT(ri — r + i-fO) J 



7=±1 



1/2 



When deriving Eq. ( |6.47 ), we used the expression similar 
to Eq. ( |6~27| ), 



(T T F(r)F(0). 



(ft(n)-ftW) ) 



(6.48) 



cos 



l(0(r- n )-0(r-T 2 )) 



All the further manipulations with Eq. (3.47) are abso- 
lutely a nalog ous to the s teps of Sec. VI A in th e deri vation 
of Eq. Q6.32| ) from Eq. fl5 28[ ). Instead of Eq. fl6.32| ), here 
we find 



G = 



dt\dt 2 



vv F E c e c J E -i (sinh7rTiisinh7rTti) 1/2 
[R A (-tr)[R A (-t 2 )r + R R (h)[R R (t 2 )}*] , (6.49) 



where the divergences (as we will see, logarithmical) 
should be cut off at times of the order of 1/Ec- No- 
tice that there is no Af dependence of the non-averaged 
conductance. The reason is similar to the absence of the 
oscillations of the capacitance in the first order in level 
spacing A - the oscillations are washed out by the quan- 
tum fluctuations of the spin mode, which is not pinned. 
We will see later that the oscillatory term in the conduc- 
tance is smaller than the leading non-oscillatory contri- 
bution to the conductance by a factor ~ (A/T) 1 / 2 . 

Let us now proceed wi th th e statistics of the elastic co- 
tunneling conductance ( 6.49] ). Ensemble averaging per- 
formed with the help of Eq. ( |6 . 2 1 a ) gives for the average 
conductance 
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(6.50) 



where C ss 0.577 is the Euler constant. At very low tem- 
peratures, T shou ld b e substituted by Aln(Ec/A), see 



discussion in Sec. VB. 

Equation J6.50| ) deserves some discussion. Firstly, we 
notice the presenc e of th e large logarithmic factor in com- 
parison with Eq. ( 6.33 ). It can be understood using the 
arguments of the orthogonality catastropheo similar to 
those applied in Ref. 27 for the inelastic cotunneling. 
Consider the elastic cotunneling process where an elec- 
tron is introduced at time t = 0, and then taken away 
at time t > Eq. Because the introduction of an elec- 
tron costs extra energy ~ Ec, all the other electrons tend 
to redistribute themselves by moving one electron charge 
through the left point contact, Fig 0. One can describe 
such a tendency as sudden change of the boundary con- 
dition in each of the spin channels. Because all the spin 
channels are symmetric, each spin mode should transfer 
charge e/2. According to the Friedel sum rule, it corre- 
sponds to the additional scattering,phase shift 5 = ±ir/2 
in each spin mode. It is knowncJ that such a sudden 
change causes creation of a large number of electron-hole 
excitations, so that the resulting state is orthogonal to 
the ground state. The probability for the system to re- 
tain its initial state during time t is P(t) ~ l/(tEc) x , 
where the index \ is related to the phase shift in all 
channels, \ = XXV* 71 ") 2 - I n our case we have four spin 
modes (two in the dot and two in the reservoir) , therefore 
X = 1. In order to find the total conductance, we have 
to sum over all possible times that electron spends in the 
dot: G cx /f-i dt/(Ect), which results in the logarithmic 



temperature dependence ( 6.50| ). The similar argument 
for the spinless electrons gives the phase shift 6 = it in 
each of the two channels. Thus, \ = 2, all the relevant 
dynamics occurs during the times smaller than 1/Ec, 
and the divergent factor disappears. 

Secondly, even though the elastic contribution is domi- 
nant in the value of the conductance, in order to find the 



temperature dependence, one h as to take into, account 
the inelastic contribution (3.44a) which yieldscilj 



G; 



Gj 



7T 3 Te- c 
' 2E C 



(6.51) 



We see that this term has stronger temperature depen- 
dence than (6.5C) and therefore, the resulting conduc- 
tance G + Gi n always grows as the temperature increases. 

Mesoscopic fluctuations of the contribution ( 6.49| ) to 
the cond uctance can be obtained with the help of 
Eq. ( PH ), 



5G(Hi)8G(H 2 



G 



-V 

2 h 



( 



In max 1 



if-. 



tit 



In(^) 



2 )V 

/ 

(6.52) 



Here the correlation field Hf is defined by Eq. ( [5.16 ), 
and the field combinations H± = H i ± H 2 are assumed 
to be much smaller than H c [the charging correlation 
field H c is given by Eq. ( p. 23 )]. The correlation func- 
tion of the conductance fluctuations starts to decrease 
fast, as 1/H±, only at fields H± > H c . Similarly to the 
case of capacitance fluctuations discussed in the previous 
Section, to obtain a representative statistics of the con- 
ductance, the magnetic field should be varied in a range 
wider than H c . 

Equation ( 6.52; ) shows that the amplitude of conduc- 
tance fluctuations are of the order of its average value, 
like in the case of the spinless fermions. However, unlike 
Eq. ( 6.34| ), the correlation function (6.52) does not reveal 
any oscillations with the gate voltage JV. 

In order to reveal this oscillatory d epend ence on Af, 
one has to expand n e ;(r ) from Eq. ( 3.44b ) up to the 
first order in action Q6.45 ). The procedure of averaging 
over all the releva nt op erators is absolutely similar to the 
derivation of Eq. (6.47), and we obtain 



ttT 



I dn . . . dT 4 i?(ri - t)R*(-t 2 )L(t 3 - r 4 ) x 
Jo 



WpEceP sin7rTr j Q 

2?t7jV ( sin7rTri sin7rT(r2 — r — 270) sin7rT(r4 — r — ijO) sin7rTr3 
sin itTt2 sin ttT(ti — t + ijQ) sin itTt4 sin 7tT(t3 — t + ijO) 



(6.53) 



7=±1 



1/2 



7T 2 T 2 



[sin7rT(Ti — t 2 + «70) sin7rT(r3 — T4 + «70) sin7rT(r3 — t 2 + 170) sin7rT(ri — T4 + «70)] 



re- 



integration over imaginary times in Eq. ( 6.53| ) is rather 
strai ghtforw ard, and tech nically very close to that in the 
Sec. |VI B 2| . Unlike Eq. fl6.49| ), here the result for the 
non-averaged conductance is A/"-dependent. We obtain 



for the oscillating contribution: 

/~i />oo 3 

r - Gr 012 / 11,11 

Uosc -vv 2 F E c eC(2n 2 Ty/ 2 J E -i ^ 



f ttT 
\ sinh irTti 



1/2-1 
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x 'ie-^R^t!) [R A (t 2 )]*L R (t 3 ) + c.c] , (6.54) 

where a 2 is a numerical coefficient: 

dxdy 



CX 2 



J sinh2/| 1 / 2 (cosha:) 3 / 2 (cosh(y — x)) 1 / 2 



11.31. 



(6.55) 



Ensemble averagi ng of Eq . ( 3.54 ) is then performed with 
the help of Eqs. (3.25) and (6.21), and the final result is 



SG osc (l)5G osc {2) 
r 2 



— as — I I In 3 ( I cos 27rn x 



A / A 



T 



(A + + A_)(A 2 + + A 2 _) 



(6.56) 



In max ( 1 ; 



A+=l- 



where the short hand notation for the conductance argu- 
ments is z = Afi,Hi, H± = Hi±H 2 , n = A/i —A/2, correla- 



tion magnetic field Hf- is given by Eq. (5.16), and the nu- 
merical coefficient 013 is given by 013 = a^e _2C /(27r 4 ) w 
0.207. 

The variance of the conductance fluctuations ( 6.56| ) in 
the unitary limit (H > H c ) is suppressed by a factor of 
four compared to its zero-field value. 



2. Finite reflection in the contact. 



For the spinless electrons, finite reflection leads to the 
oscillations in the averaged conductance already in the 
first order of perturbation theory in r -C 1. On the con- 
trary, for the electrons with spin, backscattering leads 
only to the enhancement of the oscillating part of the 
correlation function of the mesoscopic fluctuations. In 
order to demonstrate this, we expand Eq. ( 5.44bj ) up to 
the first order in the backscattering Hamiltonian (3.46). 
We obtain 



1 



E c e c ttT 
7r sin txTt 



/ dT 1 ...dr 3 .R(T 1 -T)ir(-T 2 ) V 

J ° 7=±1 



1/2 



e i2ivyM x 
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(6.57) 



sin7rTri smnT(T 2 — r — vyO) 
sin7rTr 2 sin7rT(ri - r + ijO) J [ s in7rT(Ti - t 2 + i 7 0) sin7rT(r 3 - t 2 + ryO) sin7rT(ri - t 3 + i-fO)} 



1/2' 



Performing the contour deformation for the integration 
over t 1j2 , as we did before, and the analytic continuation 
Q6.14D, we find 



G = 



Gr a 2 \r\\jEce c 
■Kvv F E c e c (2tt 2 T) 1 /2 X 

dtidt 2 T-p: > 

1- 1 (sinh7rrtisinh7rTt 2 ) 1 

[ie- i2 ^R R (ti)[R A (t 2 )Y +c.c.) 



(6.58) 



where the numeri cal co efficient a 2 is given by Eq. (3.55). 
Average of Eq. (6.58) obviously vanishes, and for the 
mesos copic fluctuations we obtain with the help of 
Eqs. (E2l 



Eq. ( 6.46] ) in the second order perturbation theory. On 
the dimensional grounds we expect this contribution to 
be 



G osc (AA) ~ G R 



\r\ 2 A 
T 



cos 2nAf. 



The low-temperature power law divergence in this f< 
mula and in Eq. ( 3 . 58) ) should be cut off at the energ; 



Ec cos 2 ttN. Calculation of the precise behavior 
of the conductance at T < e* which can be performed 
with the help of refermionization technique of Ref. is 
beyond the scope of the present paper. However, the per- 
turbation theory results indicate that the modulation of 
the conductance G osc (Af) at low temperature should be 
of the order of the average conductance G. 



5G{1)SG(2) 



M 2 A 2 
Q4 „ — m 



A ) cos 2™, 



E C T V T 

(6.59) 

where we use the same short-hand notation as in 



Eq. (6.56). The numerical coefficient 04 in Eq. (6.59) 
is given by 04 = a\e~ c /tt 4 w 0.737. 

Calculation of the contribution of the backscatter- 
ing into the average conductance requires accounting of 



VII. CONCLUSION 

In this paper, we considered mesoscopic effects in the 
Coulomb blockade (CB) regime. The emphasis was put 
on the case when the quantum dot is connected to a lead 
by a perfectly transparent single-mode channel. We have 
demonstrated that the earligc conclusion that the CB 
vanishes under this conditiont.3 is only an approximation, 
which resulted from neglecting the electron trajectories 
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returning to the channel after traversing the dot. We 
have shown that the CB persists, and its period is still 
determined by a single electron charge. However, CB 
oscillations in all the observable quantities acquire a ran- 
dom phase and therefore it is revealed in the correlation 
functions of mesoscopic fluctuations. We constructed an 
analytic, well-controlled theory describing those fluctua- 
tions. 

Our results are substantially different from the known 
results in noninteracting models of mesoscopic systems. 
For instance, the correlated ground state involves all the 
one-electron wave functions in the energy strip of the or- 
der of the charging energy Ec- The number of states in 



this strip is n ~ Ec/ 'A ^> 1 (here A is the level spacing). 
The correlation magnetic flux for the mesoscopic fluctu- 
ations <& = $o v Et I Ec is controlled by the energy scale 
Ec- The large number of states in the relevant energy 
strip leads to the robustness of the oscillatory dependence 
over about -Ec/A peaks. 

We obtained the closed analytic expression for some 
experimentally relevant characteristics. Final results are 
summarized in Tables Q for the thermodynamic quanti- 
ties (differential capacitance of an almost open dot), and 
in Table |n| for the transport quantities (tunneling con- 
ductance of an almost open dot). 



TABLE I. The ensemble-averaged differential capacitanceQ and the correlation function of the mesoscopic fluctuations in the 
Coulomb blockade regime; only the contributions oscillating with the gate voltage J\f are presented in this Table. Here |r| 2 <C 1 

1, 2 for the orthogonal and unitary ensembles respectively. For 



is the reflection coeffic ient in an al most open chann el, and (3 
more details see Eqs. (O), (O), J5.17|), and J5.22|). 





5C diff (Af)/C 


5C m {Mi)SC diff {N2)/C 2 
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3.56|r| cos 2nAf 


5.59 / A 


j cos2tt(M - Ah) 


s = 1/2 


4.53|r| 2 cos 27rAAln (^f\ 


TO) 


In (^^) + 7.12|r| 2 cos 2tt (M - A/" 2 ) 



TABLE II. Conductance of the quantum dot in a strongly asymmetric setup, see Fig. ^, and the correlation function of its 
mesoscopic fluctuations, os cillat in g wit h gate vo ltage M ■ Tun n eling condu ctanc e Gr is much smaller than e 2 /{2nK). For the 
detailed results, see Eqs. (Eggh, (|6.34|), (|6.38|), fl6.5C|) - (EH|), fl6.56|) and (Ej5§). 





G(JV)/G R 


5G(Mi)5G(M 2 )/G\ 


s = 


— (0.63 + 1.46|r| cos 2ttM) 
Ec 


0.78 / A \ 2 „ . . r . 
— [-) cos2.(AA 1 -AT 2) 


s = 1/2 


Th(?)«-^« ! )] 





Up to now, only one experiment studying the effect of 
the opening of the channel on the Coulomb blockade was 
publishedEJ. It was found that at the quantized value 
of the channel conductance, the Coulomb blockade os- 
cillations disappeared completely, in disagreement with 
our predictions. We attribute this finding to a relatively 
simple geometry of the dotcj, allowing for an adiabatic 
(rather than chaotic) propagation of an electron through 



the entire confined region. 
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APPENDIX A: DERIVATION OF EQ. ( p.7| ) 

Let us discretize the space in the direction along the 
channel axis. The fermionic Hamiltonian Hp acquires 
the form 



2ma 



(Al) 



where the transverse part of the motion is described by 
the operator 



(A2) 



and a is the discretization step. Fermionic opera- 
tors in Eq. (Al) satisfy the anticommutation relation 



{4>n(r)ipn'{r')\ = a~ 1 5 nn >5 (r± — r' ± ). The continuous 
limit of a — > 0, which will be taken in the end of the 
calculation corresponds to the usual Schrodinger equa- 
tion. Let us separate the space into two regions; re- 
gion "1" includes all the lattice sites with n < and 
region "2" includes sites with n > 0. The terms ent erin g 
into decomposition of the Hamiltonian, see also Eq. (3.6), 
Hp = H i + H 2 + H\2; have the form 



n<0 



2ma 



4>n 



Ho 



n>l 



2ma 
2ma 



(A3a) 



H 12 =- dr 



Wi + Wo 



2ma 



(A3b) 
(A3c) 



We obtain the average in Eq. (p.q) using Eq. (A3c): 



i<T T H 12 (r 1 )i7(r 2 )) 2 



(A4) 



^^{T T M^n)Mr\r 2 ))TMr\T 2 )Mr, n ). 
Am z a z 

Using the definition of the Matsubara Green function 

1 



(V>„(T)V>m(0)) 2 = Gnm(T) 



d T -H 2 



(A5) 



and low energy representation for the fermionic operator 
(only one transverse mode 4'( r ±) in the channel) 

M r -iJ = 4>o(l>{r±), 



we rewrite Eq. (A4) in the form 



~(T T H 12 {n)H(T 2 ))2 = -T t Mti)Mt 2 )x (A6) 

, 2 2 / ^^0(^)0(^)^11 (n - t 2 ; r ± , r' x ). 
4m z a z J 



As it follows from Eqs. (A3) and (A5), Green function 
satisfies the equation 



r - H j_ nn > 5 > y n > 



2ma 2 

a~ 1 5 nm 5 (ri - r 2 ) , n > 1; 



(A7a) 



2S n 



Onn' — 1 



2ma? 

a~ 1 S nm S (ri - r 2 ) , n = 1; (A7b) 



We see t hat the difference between Eq. ( A7a ) and 
Eq. ( A7b|) can be described by the boundary condition 



Gon = 0; GnO — 0. 



(A8) 



Now, we introduce coordinate x = an, and take 
the continuous limit a — > 0. Substituting G11 = 
a 2 dl x ,Q{x,x')\ x , x >=+Q into Eq. (|A6|), we obtain Eq. (E|). 



APPENDIX B: DERIVATION OF EQS. (4.14) 



We introduce more general than Eq. (4.13) correlation 
functions 

£>_ {t;xi,x 2 ) = (T T tp_(T;xi)tp-(0;x 2 )); (Bla) 
T> + (r; xt,x 2 ) = (T t ip + (t; xi)ip + (0; x 2 )); (Bib) 
D$ (t) = <T T *(r)*(0)); (Blc) 
(r; ar) - (T T $(r)^ + (0, ar)). (Bid) 



Correlation functions (4.13) are related to those from 
Eq. © by 



« , s f°° da 



da; A2?±(r; 0, x) 



x2 + x 2 "> p *+ ( T ) = p *+ (t; 0) . 



(B2) 



where the high mom enta cut-off is introduced consis- 
tently with Eqs. ( |4.l| ). Equations of motion for prop- 
agators (Bl) follow from Eqs. (4.8) and (4.10) and they 
are given by 
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(d T + iv F d x ) P_(t; x, y) = -iir5(T)sgn(x - y); (B3a) 

E c 

(d T + iv F d x )V + (r;x,y) = i sgnx — D + (r;0,y) (B3b) 

Z7T 



-■j7r<5(r)sgn(x — y); 



^2?* (r) = i^-D*+ (r;0): 

Z7T 



(B3c) 



.Ec, 



d T V^ + (t; a;) = iirS (r) - (r; 0, x) . (B3d) 

Performing imaginary time Fourier transform T>(t) = 
T En„ e" in " T r»(0„), (here ft„ = 2?mT is the bosonic 
Matsubara frequency) we find the solution of Eq. (B3b) 



V + (Q n ;x, y) = ^-{sgn(x - y) + 29 [0„(y - a:)] e^ (x ' y) } 
-^V + (n n ; 0,y) [sgnx + 29 (-Q n z) e^} . (B4) 



Substituting x = in the both sides of Eq. (B4), we find 
V + (a n ;0,y) and, then Eq. (|b|) yields: 



(fir, 



Ec 
2tt 



f 



|On|A 



(B5a) 



Propagator 2?_ is found by putting i?c* = in Eq. (B5a), 

IfiJA 



(B5b) 



and remaining propagato rs a re foun d from t he ti me 
Fourier transform of Eqs. (B3c) -(B3d) and Eq. (B5a): 



z>*(Oi) 



Ec 



2n\n n \ \n n 

it sgnf2„ 



E c 

2tt 



2?$+(f2 ra ) = 



(B5c) 
(B5d) 



In Eqs. (B5), the cut-off function f(x) is defined as 

2dy e- xy 



f{x) = 



7r 1 + y 2 



(B6) 



Inv erse Fourier transform of Eqs. (B5) gives 
Eqs. ( 4.14 ). Let us write here for completeness the result 
for the propagator 2? + (0) at arbitrary temperatures: 



T>+ (0) = In 



/ 2ttvf 
\XE c e c 



dx 



cot x 



(B7) 



At s mall t emperatures T <C Ec, Eq. (B7) reduces to 
Eq. (4.14c). At large temperatures, T ^> Ec, we obtain 



2tt 2 T 



All the res ults as sociated with the Coulomb blockade, see 
e.g. Eq. ( 4.15b| ) contain exponential terms of the form 



e ' D +(o), This gives the suppression of the charge quan- 
tization at high temperature by a factor of e _2jr T / E ° . 
This clearly indicates, that the effects considered in this 
paper and in Refs. [l4],^7| can not be obtained in any 
order of perturbation theory in charging energy Ec- 
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